A UUMPUIbH SYSTEM FOR SCHEDULING FIRE USE 



f" 1 A 
J 
/ 



I IT I I 

I H 

I . MI 



Larry S. Bradshaw and William C. Fische, 




USDA Forest Service General Technical Report INT-91 
Intermountain Forest and Range Experiment Station 
U.S. Department of Agriculture, Forest Service 



USDA Forest Service 

General Technical Report INT- 91 

February 1981 



A COMPUTER SYSTEM FOR SCHEDULING FIRE USE 

PARTI: THE SYSTEM 



Larry S. Bradshaw and William C. Fischer 



INTERMOUNTAIN FOREST AND RANGE EXPERIMENT STATION 
U.S. Department of Agriculture 

Forest Service 
Ogden, Utah 84401 



THE AUTHORS 

LARRY 5. BRADSHAW is a meteorologist for Systems for Environmental 
Management, a nonprofit research organization headquartered 
in Hissoula, Mont. He was previously employed as a mete- 
orologist by Gradient Modeling, Inc. Mr. Bradshaw earned his 
B.S. degree in meteorology from the University of Utah. 

WILLIAM C, FISCHER is a research forester for the Fire Effects and 
Use Research and Development Program at the Northern Forest 
Fire Laboratory in Missoula. His current assignment is to 
apply existing knowledge to improve fire management plans, 
with special emphasis on fire use, fuel management, and fire 
effects. Mr. Fischer earned B.S. and B.S.F. degrees from 
the University of Michigan. 



RESEARCH SUMMARY 

This report describes a user-oriented computer system that 
allows fire managers to quickly and easily analyze climatological 
data for the purpose of predicting the probable occurrence of 
desired conditions for prescribed fire. Two separate computer pro- 
grams, RXWTHR and RXBURN, make up the system. Both programs are 
designed to use National Fire Weather Data Library climatological 
data. Among the features of these programs are: adjustment of 
fuel moistures from the weather station site to the fire site; a 
newly developed duff moisture model; and the capability to simul- 
taneously consider up to 15 prescription factors. 

Step-by-step instructions for use of the programs are given 
and demonstration runs are provided. Documentation of all pro- 
grams is included. A computer terminal operator's manual for 
programs RXWTHR and RXBURN has been published as a separate Gen- 
eral Technical Report (Bradshaw and Fischer 1981) . 
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INTRODUCTION 

The timely occurrence of desired burning conditions has been and continues to be a major 
source of uncertainty associated with using fire as a resource management tool. Although such 
uncertainty cannot be eliminated, it can be reduced by using information contained in clima- 
tological records (Fischer 1978); for example, analysis of climatological data can describe an 
area's average weather and also likely deviations from this average situation, A fire manager 
can use such information to identify periods when desired burning conditions might occur. 
Further analysis can result in calculation of the probability of desired fire conditions 
occurring during the periods identified (Furman 1979). 

The value of cliniatological analysis as an aid for scheduling prescribed fires can be 
summarized as follows: 

1. Prescribed 'fire schedules can be geared to the number of probable burning days. Over- 
schecluling of fires, especially during traditional burning periods, can be minimized. 

2. Prescribed fire opportunities outside traditional burning periods can be identified 
insofar as available weather data will allow. 

3. Project priorities can be more nearly matched to the probable occurrence of prescribed 
conditions. Some fires require a narrow range of burning conditions for successful accomplish- 
ment. These fires can be given top priority for scheduling during periods when such conditions 
are most likely to occur. 

4. Fire management prescriptions for preplanned use of chance wildfires can be evaluated 
before they are implemented. The frequency and time of occurrence of prescribed fire conditions 
can be identified. This information can then be used to evaluate the probable impact of the 
prescription and then compared with the probable performance of alternative prescriptions. 

Despite its obvious value, fire managers rarely have had the means to easily analyse 
climatological data when planning fire use. This is due, in part at least, to the tedious 
nature of the task. It is also due, in part, to weather data problems: missing data, lack 
of data during spring and fall, and, in some cases, unreliable data resulting from poor 
observation and poor weather station maintenance. 

This report describes a user-oriented computer system that allows fire managers to quickly 
and easily analyze climatological data for the purpose of predicting the probable occurrence of 
desired prescribed fire conditions. Instructions are provided to assist both the fire manager 
and the computer specialist through the various steps required to use the system. Finally, 
documentation of the computer programs used is provided. A computer terminal operator's manual 
has been published as a separate report (Bradshaw and Fischer 1981). 

SYSTEM DESCRIPTION 

Two separate, but interrelated computer programs, RXWTIIR (Prescribed Fire Weather) and 
RXBURN (Prescribed Fire Conditions) make up the system. Both programs are designed to use the 
climatological data base stored in the National Fire Weather Data Library (NFWDL) (Furman and 
Brink 1975), which is located at the USDA Forest Service's Fort Collins Computer Center (FCCC) , 
RXWTIIR and RXBURN are stored in the Region 1 program library at FCCC and are available via remote 
terminal or batch processing to all who have access to FCCC facilities. 

Information Provided by the programs 

Program RXWTHR provides climatological summaries and co-occurrence frequencies of user- 
selected fire weather and fire danger rating parameters. The 16 parameters available for 
analysis are: 



1. State of the weather 

2. Temperature 

3. Relative humidity 

4. Wind direction 

5. Wind speed 

6. Maximum temperature (last 24 hours) 

7. Minimum temperature (last 24 hours) 

8. Maximum relative humidity (last 24 hours) 

9. Minimum relative humidity (last 24 hours) 

10. Precipitation duration 

11. Precipitation amount (last 24 hours) 

12. 1 hour fuel moisture 

13. 10 hour fuel moisture 

14. 197S National Tire Danger Rating System (NFDRS) Energy Release Component (ERC) 
(with selected fuel models) 

15. NFDRS Burning Index (BI) (with selected fuel models) 

16. Duff Moisture (24-hour average for entire layer) 

The user is reminded that Program RXWTHR can be used to summarize and analyze only those 
parameters that are included in NFWDL for the weather station of interest. The importance of 
maintaining complete records at the NFWDL and keeping them current through the Adminstrative 
and Forest Fire Information Retrieval and Management System (AFFIRMS) (Helfman and others 1975) 
cannot be overly stressed. 

Examples of the type of information provided by RXTHWR are shown in figures 1, 2, 3, and 4, 
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Figure 1. --Example of Program RXWTHR summary table for temperature. 
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Figure 2. --Example of Program RXWTHR summary table for wind direction, 
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Figure 3. --Example of Program RXWTHR two-way co-occurrence table for wind direction and 

wind speed . 
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Figure 4 . --Example of Program RXWTI1R three-way co-occurrence table for temperature, relative 

humidity, and windspeed. 

Program RXBURN provides analyses of local prescription condition frequencies based on 
user-defined prescription conditions. The 16 parameters available for use as prescription 
conditions are : 

1. State of the weather 

2. Temperature 

3. Relative humidity 

4. Wind direction 



5 , Wind speed 



0. Maximum temperature (last 24 hours) 

7. Minimum temperature (last 24 hours) 

S. MaAimum relative humidity (last 24 hours) 

9. Minimum relative humidity (last 24 hours) 

10. .Number of days since last measurable precipitation 

11. Precipitation amount (last 24 hours) 

12. 1 hour fuel moisture 

13. 10 hour fuel moisture 

14. 1978 NFDRS ERG (with selected NFDRS fuel models) 

15. 1978 NFDRS BI (with selected NFDRS fuel models) 

16. Duff moisture (24-hour average for entire layer). 

As many as 15 of these parameters may be used in a single prescription. Information pro- 
vided by RXBURN is exemplified in figures 5, 6, 7, and 8. 
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Figure 5. -Example of Program RXBURN summary table, 
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PROGRAM RXBU7NI RUN MO. 1 DISTRICT! FIRELA9 



FOREST! MONTANA 



PAGE NO. 2 



PRESCRIPTION OCCURRENCE RY 10 DAY PERIOD AMD MONTH 




* MONTH 


PERIOD 
BEGINS 


PREFERABLE DAYS ACCEPTABLE DAYS UNACCEPTABLE OAYS * 
NO. DAYS MEAN NUMBER ERCENT MiAN NUMBER PERCENT MEAN NUMBER PERCENT 


* MAY 


1 







5 


5K 


1 


18 


17* 


5 


81 


7fl% 


* MAY 


11 


109 





5 


5* 


2 


33 


30% 


4 


71 


65% 


* MAY 


21 


118 


1 


8 


7* 


3 


40 


34% 


4 


70 


59% * 


*> MAY 


TOTAL 


331 


1 


18 


5* 


6 


91 


27* 


14 


222 


67% * 


* 

























JUN 


1 


110 


t ft tt toft ft i 

l 


11 


#ft< 

10* 


3 


t ft*4HH 

50 


45* 


3 


49 


45% * 


JUN 


11 


109 





4 


4JS 


2 


39 


36* 


4 


66 


61% * 


* JUN 


21 


106 





6 


6% 


3 


45 


42* 


3 


55 


52% * 


* JUN 


TOTAL 


32 5 


1 


21 


6* 


a 


134 


41* 


11 


170 


52% * 


* JUL 


1 


13 9 





#* 
3 


tut** 
2% 


5 


78 


*! 

57* 


HtflB 

4 


<tttttt 
57 


><*** 
41* 


JUL 


11 


150 


1 


8 


5* 


4 


62 


41* 


5 


80 


53% 


* JUL 


21 


165 





5 


3* 


4 


63 


38% 


6 


97 


59% 


* JUL 


TOTAL 


453 


1 


16 


4% 


13 


203 


45% 


15 


234 


52% 




























##***ffft#ft#ft*nnntft##ftft#**#* 


*##*ftft#*ft* # *## * ft* * *ttft ft ft * 


* AUG 


1 


159 





7 


4* 


4 


66 


42% 


5 


06 


54% * 


* AUG 


11 


146 





3 


2* 


4 


56 


38% 


5 


87 


60% 


AUS 


21 


16o 


1 


10 


6K 


4 


67 


42* 


5 


83 


52% 


AUS 


TOTAL 


465 


1 


20 


4% 


12 


189 


41% 


16 


256 


55% 


* 

























#tt<ttt*fftt9***t***#t***#*B*e##tttt#**#BO*tl*Oti 


SE 3 


1 


137 


1 


12 


9JS 


3 


55 


40% 


4 


70 


51% 


SEP 


11 


129 


1 


11 


9 


3 


45 


35% 


5 


73 


57% 


* SEP 


21 


117 


1 


11 


9* 


2 


39 


33* 


4 


67 


57% 


I 5ED 


TOTAL 


383 


2 


34 


9* 


9 


139 


36% 


13 


210 


5S% 


#< 

OCT 


ntt*ft**tt 
1 


***#( 
90 


1 


t*tt*fttt* 

9 


#< 
10% 


HHHtttHHtl 

1 


HJBtt##l 

20 


t*tt 
22% 


4 


ftfttttte 
61 


lttO*lt 

68% 


* OCT 


11 


90 





3 


3* 


1 


21 


23% 


4 


66 


73* 


OCT 


21 


78 








0* 


1 


22 


28* 


4 


56 


72% 


OCT 


TOTAL 


258 


1 


12 


5* 


4 


63 


24% 


11 


183 


71% 


* 

























* it*** ft* fti 
NOV 


l U * &* 

TOTAL 


*** 



!** 




h*HllH***(KH* 




* ** 
0% 


Ht(HHHHtt 




tft*Htttttft 




hi-ttfttt44ttf 

0* 


nm *** 



& *imin*flfl 



^^l tt^-tttf^ttti 

OB 

n 


*#,,),#*), ttt*ftOftftftftttftft*ftft*ftftft*tt 


TOTAL 


TOTAL 


2215 


8 


121 


5* 


51 


819 


37* 


80 


1275 


58% * 


* 






















* 



Figure 6. Example of Program RXBURH prescription frequency of occurrence table 



RXBlPNi *UM -JO. 1 DISTRICT! FIRELAB 

PRESCRIPTION RUN LF^GTH 



PAGE 



PREFERABLE DAY RUNS * ACCE"Tft3LF DAY RUNS 
PERtOD * OERCENTKES * "ERCFNTILES 

HEANJ 25TH MEDIAN 75TH MEA^J ?5TH MFDIAN 7iTH 



J^ACCEPTARLE OAY RUNS 

PEHC^MTILES 
MEAN ?5TH MF31AN 7HTH 



* -..._ 


1 


3 


1 


1 


4 1 


1 


I 







2~ 


----- 


7 


* HAf 


11 


1 


1 


1 


1 ? 


1 


2 


? 


3 


1 


P 


S 


*Af 


21 


1 


1 


1 


1 ? 


1 


1 


? 


3 


1 


? 


5 


; 


TOTAL 


1 


1 


1 


1 ? 1 1 

)(HHKt&0-|Hftt*HtflOttflltttfttt*ttttUtttt& 


tf tt ft fi H 


* 


1 


3 


ft 


* JUN 


1 


2 


1 


1 


2 ? 1 1 


3 


* tt # tt 1 

? 


l H (t 
1 


IB rtitt 1 
? 


3 


* J'JY 


11 


1 


1 


1 


1 ? 


2 


? 


3 


1 


1 


? 




JJV 


21 


1 


1 


1 


1 ? 


1 


1 


2 


3 


1 





, t 


JJ>J 


TOTAL 


1 


1 


1 


\ 5 


1 


2 


2 


3 


1 


? 


> 



* 


JJ. 


1 1 


1 1 1 


? 1 


uttnnn 
2 


3 


o it H ff i 
2 1 


1 P ft 


i 
* 


JUL 


n t 

21 1 


1 1 1 
1 1 1 


2 1 


J 


2 


2 1 
1 i 


* r " 

3 


* 


JJL 


TOTAL 1 


1 1 1 


? 1 


1 


2 


j i 
2 1 


? ft 
2 3 




* 


ALJ5 


1 1 


1 1 1 


o # # a # * u * 


UttUtt* 


ntinnnnnt 
p 


***# 


" *# 


* 


AU5 


11 1 


1 1 1 


? 1 


1 


L. 

3 


3 1 


? 3 * 


* 


AU3 


21 1 


1 1 1 


? \ 


J 


? 


"* i 

? 1 


i 


t 


AU3 


TOTAL 1 


1 1 1 


? \ 


1 


) 


^ I 
1 1 


2 3 


i 














i 1 


2 3 

n 










OlfHI(ltMttlltt tt ttVttttHn 




* 


s 


1 1 
11 1 
21 1 
TOTAL 1 


1 1 1 
1 I 1 

1 1 1 
1 1 1 


? 1 

? 1 
? 1 
? 1 


1 

2 
1 
1 


2 
1 
? 
2 


2 1 
3 1 
3 1 
3 1 


2 2 a 

y <3 * 

3 3 
J? A <t 


i 


OCT 


1 1 


1 1 i 


'#*tt*IHHHH* # 
11'- 



fttfttttftttttttittft 


* 
* 


QCT 
3CT 

ocr 


11 2 
2\ 

TOTAL 1 


1 1 2 



1 1 1 


t 1 

? 1 
? 1 

3 1 


j 
1 
1 


3 
? 


* 2 
2 1 


3 6 

1 ? 










C 1 


1 


? 


* 1 


? 5 






-JOV 


TOTAL 


000 


# * tt t ttoit 



mntB 



IB * K H 

OOQ 


H 
H HttHHtHllll 





7,-Ex^le of Progra, RXBURN proscription run loilB th sundry. 



PROGRAM RxaU'N! W4 W. 1 DISTRICT: FIRF.I.AR r <HEST: MONTANA 

PROB43ILITY Or MFETING PRESCRIPTION 1. 2, AVD 3 DAYS 14 THE FUTURE 



PAGE NO 







MONTH: 


MAY 








& 




















TODAYS 


TOMORROW ,* 


# 


2 DAY 1 ! 


o** 


** 3 DAYS 


* M DAYS 


* 


COND 


PREF ftCCP UNAC 


PREF 


ACCP 


UNAC 


PRFf fiCCP 


JNAC 


o 
















D 


3REF 


17* 50* 33* 


9% 


38* 


S3fc 


7% 32* 


61* PREF S 


B 


ACCP 


12* 47* 4U 


8% 


35* 


S&* 


7<6 31% 


62* ACCP 85 





JNAC 


2* 18* 80* 


4* 


24, 


72* 


5% 26* 


6B* UNAC 216 


B 







































MONTH! 


JUN 




























TODAYS 


* TOMORROW ** 


## 


2 DAYS 


* 


*** 3 DAYS 


* ** DAYS 


tt 


COsjO 


PREF ACCP UNAC 


PREF 


ACCP 


UNAC 


BRFF ACCP 


JNAC 



















a 


3REF 


20* 60* 20* 


10% 


SI* 


39* 


B% 45% 


47* P=EF 20 


B 
















H 


ACCP 


9* 56* 35% 


8* 


47* 


4i% 


7* 4355 


49* ACCP 132 


* 


JNAC 


3S ?8* 69* 


54 


37* 


53* 


6% 40% 


5418 tmc 173 


H 
















1 
B- ji n 
















T W If 






MONTH: 


JUL 


tt 19 frtt ff ft 






i> 


TODAYS 


TOMORROW ** 


* 


2 DAYS 


B 


" 3 DAYS 


* n DAYS 





CONO 


PREF ACCP UNAC 


PREF 


ACCP 


UNAC 


PREF ACCP 


JMAC 





'REF 


0* 75* 25* 


4* 


48* 


4^* 


4 45% 


51* PREF 16 


ft 


ACCP 


4* 52* 44* 


4* 


46* 


51% 


4* 45% 


52* ACCP 206 





JNAC 


3* 36* 61* 


3* 


43* 


54* 


<<& 44% 


52* UNAC 227 


ft 



Figure 8. --Example of program RXBUUN prescription occurrence contingency table 



Information Needed to Use the Programs 

Both programs require climatological data from the National Fire Weather Data Library. 
Information on the NFWDL data base structure and use is detailed by Furman and Brink (1975), 
For completeness, instructions for accessing NFWDL data for use in RXWTHR and RXBURN are outlined 
in the User's Guide in the next section of this report. 

The user must also furnish information requested on the information worksheets for RXWTHR 
and RXBURN (figs. 9 and 10). 



r i re >>eather : 


it at ion [nfonut ion 


Stat ion name 


, N'O. 


Elevation 


ft .Latitude 


.Climate class--' 


.Slope class 


Fuel model- 


.Last frost 


- .Grass type ; annual 


, perennial 


"i ear begin 


, 'lear end 


.Date beiun 


.Date end 


SLZ Ad-u^-aer 


it Factors i if anvl 


= south, 4=west ], Site elevation 


r t . 


Asrect i 


l=north, >east, 5 


Canopv cover 


1 L=open, 


>elosed ) 





Horizon Informatun if Duff Moisture selected. See User's Guide, aj>pj?ndi.x 0) : 



Duf^joil Type Thickness 



4 



cm 
cm 
cm 
cm 



cm 



\J Use up to 80 characters 

2J See User's Guide, appendix C 

3/ See User's Guide, appendix B 

4_/ For NFDRS indices only. See User's Guide, 
appendix A 

5/ For NFDRS indices only 

* + *T* ******-******#***** ff + **^ Al>t .^ # ^ 



Summary Tableisj Requested (Select up to 15] : 

State of the weather 

Temperature (degrees F) 

Relative humidity (0 

'-und direction (S point) 

lund speed (mi/h) 



_Ma.< temperature C4 h, deg. F) 



_Min cemnerature i 14 h, deg. F) 
_Ma\ relative humiditv (21 h, J ,l 



_Min relative humidity (.24 h, "i 
_Precip duration (lust 24 h) 
_Precip amount (24 h, 0.01 in) 
_1 hour fuel moisture ( D o) 
_10 hour fuel moisture ("a) 
_NFDRS ERG 
_\'FDRS BI 
_tXit"f Moisture (24 h average, *, 



jO-occurrence Table) si Requested (If ''-uav 
selected, Wind Direction 


table desired leave last space blank. If 


must always be isted fir<*t ") 


with 


with 


with 




^^-^ ^^^^ 

iv 1 1 h 


. with 


- ...-. 

with 


. with 


'.vith 


with 





Figure 9.-RX1VTHR user information sheet. 
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Total lumber Different Stations in This Run 



User's Name .Subunit .Unit 

Project^ 

Fire Weather Stat ion Informat ion : Station name 

Elevation ft . , Latitude .Climate class^ .Slope class- 



Fuel model .Last frost- .Grass type^- ; annual .perennial 

\ear begin . .'tear end .Date bevjin _,_Date end 



Site Adjustment Factors (if any)- 

Aspect ( l=north, 2 = east, 3-south, 4=west ), Site elevation ft . 

Canopy cover ( l=open, 2 = closed ) 

Duff/Soil Horizon Information (if Duff Moisture selected. See User's Guide, \ppendix D) : 
Layer Duff/Soil Type Thickness 

1 cm 

2 cm 



5 cm 

4 cm 



cm 



\J Use up to 80 characters 

2f See User's Guide, appendix C 

3/ See User's Guide, appendix B 

4_/ For MFD PS indicies only. See User's Guide, 

appendix A 
5/ For NFDRS indices only 



****************** + ****** 



Prescription Factor Selections (Check and set limits for up to 15 factors. 1: 

Preferable Rx Acceptable Rx 
Factor Minimum Maximum Minimum Max i mi 



State of the weather 
Temperature (deg . F) 
Relative humidity (V 



Wind direction (3 point) 

Wind speed (mi/h) 

Max temperature (24 h, deg . F) , 
Mm temperature (.24 h, deg. F) . 
Max relative humidity (24 h, 
Min relative humidity (.24 h, ) 
Days since last precipitation.. 
Precip amount (24 h, 0.01 in).. 

1 hour fuel moisture (%) 

10 hour fuel moisture (%) 



_NFDRS ERG 

_NFDRS BI 

Duff moisture (24 h average, "0 



Figure 10.--RXBURN user information sheet. 



Program Use 

Ideally, the programs should be used sequentially. RXWTHR will give the manager a method 
for screening possible combinations of prescription factors to eliminate those which have a 
low probability of occurring simultaneously, as well as supplying a local fire climatological 
summary, 

Once a feeling for the occurrence and timing of desired prescription conditions is reached, 
RXBURN may be used to provide a detailed summary of prescription occurrence frequencies. 

Because of the costs involved, RXWTHR should be run as few times as possible for any given 
station, while RXBURN may be run many times for different prescriptions and planning processes. 
Each program is, however, completely independent of the other and sequential use is not 
mandatory. 

Options and Features of the Programs 

FORTRAN IV code is ASCII standard, 

1978 National Fire Danger Rating System (NFDRS) algorithms and subroutines are used by 
the programs where applicable (Deeming and others 1977). 

The routines operate on simple input streams. 

The programs contain input error checks that print error messages. 

The programs offer 1978 NFDRS indices, Energy Release Component (ERC) and Burning Index 
(BI) as prescription parameters based on NFDRS fuel models. 

The programs contain site adjustment factors to adjust fuel moisture values (and 
resulting ERC's and El's) to locations with elevations, aspects, and canopy covers different 
than those of the user- specified station. 

RXhTHR allows as many as 15 summary tables and five co-occurrence tables to be produced 
in a single run. 

RXBURN allows as many as 15 prescription factors to be simultaneously considered in a 
prescription condition, 

As many as 99 fire weather stations may be analyzed for summaries or prescription 
occurrence in a single run. 

Repeated runs of the same fire weather station allowing different prescription con- 
ditions to be analyzed in one run is possible through the construction of a disk file to hold 
rirst nm computations. 

A newly developed duff moisture model is available (Fosberg 1975). 



factor One U is l^TJ ^T^T, Condition ranges to be entered for each prescription 
burning condition " preferable " bu in * condition, the other is an "acceptable" 

RXBURN allows the number of days since last precipitation as a prescription condition. 
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System Limitations 



Programs RXWTHR and RXBUR.N are designed to use archival data in the National Fire Weather 
Data Library (Fur man and Brink 1975). Gibson (1977) and others have noted a limitation in the 
use of National Fire Weather Data Library records for planning prescribed f ireobservations are 
taken only once a day over the length of the fire season. This results in two constraints on any 
system that uses the data. First > the single observation must be considered to be representative 
of the entire day, when in fact it represents the "typical worst" conditions (midafternoon,, 
southwest aspect, midslope, open canopy). Second, weather observations are rarely recorded for 
more than 5 months each year during the "fire season." Data therefore tend to be scarce in 
prefire and postfire season periodstimes that prescribed fire operations are generally most 
active. 

To reduce computer storage requirements, program RXWTHR analyses a maximum of 5 months of 
weather data from the National Fire Weather Data Library in a single run. This is usually not a 
limitation of the system since few stations have more than 5 months of data on file for each 
year. For those stations with more than 5 months of fire weather data, two or three runs of 
RXWTHR may be needed to obtain the desired climatological summaries, The 5 -month restriction 
does not apply to program RXBURN. 

A newly developed duff moisture model is available for test use (Fosberg 1975). The model 
is based on theoretical considerations of water transport through the litter, duff, and soil 
horizons as affected by daily weather observations. The model is still in the validation stage, 
and results should be viewed as tentative and used with caution. Preliminary results indicate 
the model will consistently underestimate duff moisture content, 

USER'S GUIDE 

Step-by-step instructions for using RXWTHR and RXBURN follow. Demonstration runs of both 
programs are included to support these instructions. The instructions are written so that anyone 
familiar with the use of remote terminal or batch processing equipment should be able to success- 
fully run the programs, 

In organizations that employ computer technicians, the fire manager's role Ln running these 
programs may be limited to completing the user information sheets (figs. 9 and 10). 

Create a Data File 

Two items of information are needed to obtain data from the NFWDL. First is the six-digit 
code (or codes) of the fire weather station(s) to be analysed and the years of data to be ana- 
lyzed. Second is the file name in the library that contains the lowest station code that will be 
used in the analysis. For example, if the stations to be analyzed are 034567, 245789, and 
003452, only the file name that contains station 003452 is needed. 

Obtain a File Name 

There is a possibility that a local fire weather officer or computer specialist has a 
current listing of the files in the library. If not, obtain a file list in the following manner; 

12345 

12345678901234567890123456789012345678901234567890 



@ASG,A FIREDATALIB*PROGRAMS. 

FIREDATALIB*PROGRAMS.LISTFILES 
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System software will then respond with a listing of file names and stations in the file. An 
le of the general format follows. (Note: ssssss represents the six-digit station code, and 

vy represents the last two digits of the year that data begins (FROM) or ends (THROUGH); nn, rain, 

and oo represent assorted numbers and letters of file names.) 

FILE NAME STATION YEAR LIMITS DATE OF LAST 

FROM THROUGH UPDATE 

riREDATALIB*nn-mm ssssssyy ssssssyy mmddyy 

FiREDATALIB*oo ssssssyy ssssssyy mmddyy 

Scan the station year limits column until the group containing the lowest six-digit station 
code of stations to be analyzed is found. The entire file name (FIREDATALIB*mm-mi) is to be used 
in place of "FILE" in the following data acquisition sequence. 

Create a Card Image File 

In creating a card image file for use in RXWTHR or RXBURN, it is wise to generate a user 
program file of the data. This allows the data file to be stored in the Mass File Directory at 
fCCC for 6 days from the date of creation, allowing subsequent runs of RXWTI1R or RXBURN to access 
the data with recreation of the file on each run. This is particularly helpful in the event of 
input stream errors resulting in program termination prior to complete execution, 

This process is accomplished by executing the following control sequence at the PCCC . 
(Note; <\gain, ssssss represents the six-digit station code; yy the year to begin data inclusion 
in the file; and the second yy the year to end data inclusion. If all available years are 
requested, use yy = 00 and yy = 99.) 

1234 
1234567890123456789012345678901234567890 

3RUN,... 

QASG,A FIREDATALIB*PROGRAMS. 
@ASG,A FILE. 
QUSE 2,, FILE. 

QASG.UT NAME. (NAME MAY BE ANY NAME) 

@USE 15. , NAME. (USE SAME NAME AS PREVIOUS LINE) 

FIREDATALIB*PROGRAMS.GETDATA2 

ssssssyy ssssssyy 

etc., until all stations are listed in ascending order. 



FREE 2. 

3ASG,T 2. < THIS ASSIGNS TEMPORARY DISK FILE USED IN PROGRAMS) 



o H 4" f the Statl n numbers d the nu.ber of card 

-cad wea; da f' loca unit 15 " an3lySiS ^ ^^ RX R r RXBURN ' which 
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Subsequent runs using the same station within the next 6 days would use the following 
sequence in preparing the data file for program use: 

1 2 

12345678901234567890 

@RUN,. . . 

@ASG,A NAME. (NAME IS ONE PROM ORIGINAL RUN) 

@USE 15. ,NAME. 

@ASG,T 2. 

Using Program RXWTHR 

To use RXWTHR, four steps must be followed. 
Step 1. Fill out user information sheet(s). (See fig, 9,3 

One sheet should be completed for each station to be summarized, or for each summary for 
a station. The entries on the sheet are mostly self-explanatory, and are briefly summarized 
here for completeness. 

a. Total number of stations in run. --This value is the number of different stations that 
will have summaries produced in one run of RXWTHR. 

b. Fuel model, grass type, and_date ojf last spring frost . --These values need only be 
entered when NFDRS indices are requested in output tables, NFDRS fuel model descriptions are 
detailed in appendix A, 

Ct Slope class . --This value must be entered if NFDRS indices or site adjustment options 
are selected. Slope classes are defined in appendix B, 

d. NFDRS cl imat_e___class . - -This value should be entered on all information sheets and as 
detailed in appendix C. 

o. Activity information. --User comments and documentation are entered here. As many as 
80 characters are allowed. Suggested entries include user name, project objectives, site 
adjustment factors, district, and other information. 

f. Duff/ soil horizon information. --This information is only required by the selection 
of duff moisture as an output parameter. Instructions are detailed in appendix D. 

g. Summary _. t_abl_es_. --If summary tables are desired, select parameters for inclusion. 

h, Co -occur r once tables. --If co-occurrence tables are desired, select up to five by 
filling the blanks with the desired parameters from the total list in the preceding summary 
section. The only limitation is that wind direction, if selected, must be the first param- 
eter entered in a table composition. 

Step 2. Create data file for program use and load program RXWTHR, 

File creation was previously covered; program loading is done by the execution of the 
following control cards: 

123 

123456789012345678901234567890 

@ASG,A CSSG*R1LIB. 

@XQT CSSG*R1LIB. RXWTHR 
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Step 3. Transfer information from user sheet(s) to machine readable formats. 

The information from each information sheet constitutes a directive block that tells the 
program what options are requested and supplies other needed information. There must be one 
directive block for each station in the analysis and for each multiple run of a single station. 
Information entered in one directive block does not need repeating in subsequent station analysis 
nor multiple runs if the information does not change. Details of options and information entry are 
detailed below with examples. 

a. Total number of stations card. --This card must be the first in each directive block. 
The value" is entered in columns 1 and 2, right justified, no decimal. 



1234567890 



nn 



b. S t at i o n , years , and dates of an a 1 y s i s cards . - - Th e first card in this sequence instructs 
the program that three cards of station, year, and dates information follow. It consists of 
the word STATION beginning in the first column. The first trailing card contains complete 
station information as detailed in table 1. This is then followed by two more trailing cards 
with additional information, The first identifies the years to begin and end data inclusion 
to the analysis, the second identifies the seasonal dates of data inclusion. If all available 
years are desired, use 1900 to 1999. Program RXWTHR is restricted to 5 months of data analysis 
in a single run. For stations with more than 5 months of archival data, see the following 
section on System Options. For an example, consider Philipsburg Ranger Station for the sea- 
sonal dates of June 1 to September 15 for 1960 to 1977. 



1 23456 7 

1234567890123456789012345678901234567890123456789012345678901234567890 

STATION 

PHILIPSBURG RS 243002 5280 46 H 3 P 3 0615 T 3 6200 1 

YEARS I960 1977 

DATES 0601 0915 

c. Activity information^card..- -The first card in this sequence tells the program that 
one card of activity information follows. It consists of the word ACTIVITY beginning in 
column 1. It is followed by a card that contains the information under activity on the user 
information sheet. 



123456 7 8 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

ACTIVITY 

FIRE WEATHER SUMMARIES FOR PHILIPSBURG, LARRY JONES, PROJECT 5 
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Table 1. --Station information card 



Input variable 


Limits of 
value ranges 


Card columns 
(inclusive) 


Format 


Information 
source 


Station name 


15 characters 


1-15 


4A4 


user 


Station number (6 digits) 


0-999999 


17-22 


16 


user 


Station elevation (ft) 


0-99999 


24-28 


F5.0 


user 


Station latitude (degrees) 


0-99 


30-31 


12 


user 


NFDRS fuel model 


A-U (except M) 


33 


Al 


appendix A 


NFDRS slope class 


1-5 


35 


11 


appendix B 


Grass type 


A or P 


37 


Al 


appendix A 


NFDRS climate class 


1-4 


39 


11 


appendix C 


Date of last frost 


0101-1231 


41-44 


14 


user 


Repeat run ? (REPEAT) 1 


T or blank 


46 


LI 


user 


Save run ? (SAVE) 1 


T or blank 


48 


LI 


user 


Site adjustments ? (ADJUST) 1 


T or blank 


50 


LI 


user 


Site aspect code 


1-4 


52 


11 


user 


Site elevation (ft) 


0-99999 


54-58 


F5.0 


user 


Site canopy cover 


1-2 


60 


11 


user 



If the run is to save computations for the next run to use, enter a T in card column 48, other- 
wise leave column 48 blank. 

If the run is to use computations from the last run, enter a T in card column 46, otherwise 
leave column 46 blank. 

It the run is to make site adjustments, enter a T in card column 50. Be sure if site adjusments 

are to be made that correct values are entered for aspect, elevation, canopy, and NFDRS slope 

class. 

If the run is not to mak'e site adjustments, leave columns 50 through 60 blank, 

logical variables. 

d. Duff/soil information. --This input sequence begins with a card instructing the pro- 
gram as to how many layers (cards) of duff/soil information follow. It consists of the word 
DUFF in columns 1 through 4, and the total number of layers (2 to 5) in column 12. No decimal 
is punched. This card is then followed by as many as five cards that define the layer type 
and thickness (in centimeters) . Layer type is in columns 1 through 5 and thickness in columns 
6 through 10. Both values have decimals punched. 



1 2 

12345678901234567890 

DUFF 3 

1. 3.5 

2. 7.0 
7. 300.0 
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e. Summary table option cards. --This sequence begins with a card that instructs the 
program that summary table options follow. The word SUMMARY is printed beginning in column 
1 and is followed by as many as IS cards defining the parameters to be summarized. There is 
one card for each parameter checked on the sheet. The entire parameter name should be spelled 
correctly beginning in the sixth column of the card. Information in the parentheses may be 
excluded. After the last option card has been printed, instruct the program that the last 
card has been read by placing the word END beginning in the first column of the next card. 



123 

123456789012345678901234567890 

SUMMARY 

TEMPERATURE 

WIND SPEED 
END 



f. Co-occurre-nce tab! e jjption cards . - -Thi s sequence begins by instructing the program 
that co-occurrence option cards follow. This is done by use of the word CO-OCCUR beginning 
in the first column. It is followed by as many as five option cards that define what tables are 
to be produced. Values are transferred from the user sheet to cards, one card per table. 
Make sure the spelling is the same as the parameters in the summary list. The first param- 
eter begins in the sixth column., the second in the 31st, and third (if any) begins in the 
56th column. Again, information in the parentheses should be excluded and the last option 
card is followed by an END card. 



1234 567 
1234567890123456789012345678901234567890123456789012345678901234567890 



CO-OCCUR 

WIND DIRECTION 
TEMPERATURE 

END 



WIND SPEED 
RELATIVE HUMIDITY 



WIND SPEED 



g. The RUN card. --This card must be the last card in each directive block. It consists 
simply of the word RUN in columns 1 through 3 and signals that all cards in the block have been 
read. The program then begins execution. 



1 
1234567890 

RUN 



Step 4, Terminate program execution. 

The program will continue execution, reading one directive block at a time, and perform- 
ing the requested calculations until coming upon an end-of-file mark (@EOF) when reading the 
first card of the next directive block (Number of Stations). After all directive blocks have 
been entered, the program is terminated by the following commands; 
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1 

1234567890 

@EOF 

@FIN 



Program RXWTHR Error Messages 

1. Control sequence incorrect, program aborted. 

Cause; A control card was read that was uninterpretable by the program's internal diction- 
ary. Most common cause of an aborted program is misspelled control words or not 
starting the word in column 1 of the control card. Can also be caused by not 
placing an END card at the end of an option series. 

2. No summary or co-occurrence option specified, program aborted. 

Cause: There is neither a SUMMARY nor a CO-OCCUR control card in the input stream. 

3. Station numbers on descriptor card and data file do not match, program aborted. 

Cause: Keypunch error on STATION card following STATION control card, or error in data 
file structure. Remember that stations must be analyzed in ascending order. 

4. At least one summary option card is incorrect, program aborted. 

Cause: One of the option cards trailing the SUMMARY control card is misspelled or does 
not begin in column 6. 

5. At least one co-occurrence option card is incorrect, program aborted. 

Cause: One of the option cards trailing the CO-OCCUR control card has a misspelled word, 
or the words do not begin in the correct columns. The correct columns are 6 and 
31 for two-way tables, and 6, 31, and 56 for three-way tables. 

6. Difference between site and AFFIRMS station greater than 2,000 feet, no adjustments made. 

Cause: When using ADJUST option, elevation differences are limited to 2,000 feet. Program 
executes, but without any site adjustments. 

7. More than 5 months weather summary requested, only first 5 months processed. 
Cause: Self-explanatorysplit data into two or three sections. 

Program RXWTHR Output Quantity 

For each station, the following output is produced: 

SUMMARY tables: One (1) page for each parameter requested. 
CO-OCCURRENCE tables: 

Two-way: One (1) page for each table per 2 months analysis; for example, a two-way 
table from May through August will have a two-page output. 

Three-way: One (1) page for each table for each month's .analysis; for example, a three- 
way table from May through September will produce five pages of output. 
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Run priority 




Cost 






CAU 


(sec) 


DEMAND 


$6. 


00 - 


$7. 


00 


14 - 


15 


L 


S3. 


00 - 


$4. 


00 


14 - 


15 


M 


52. 


00 - 


$3. 


00 


14 - 


15 


P 


$1. 


00 - 


$2. 


00 


14 - 


15 



Program RXWTHR Cost Estimates 

lor A tvpical RXKTIIR run (one station, 5-month analysis, five summary tables, two to 
three co-occurrence tables) a user can expect costs and resource use as outlined below; 

Resource time (sec) 
22 - 25 
22 - 25 
22 - 25 
22 - 25 

RXWTHR Demonstration 

Putting together steps 1 through 4 results in the execution of a very simple run of 
RXI'.TIIR, for demonstration purposes, a sample analysis of two stations in a single run will 
be used. Other runs are demonstrated in the Program Options section. In tins example, 
Philipsburg Ranger Station and West Glacier headquarters are the two stations in the sample. 
Assume data file is FIREDATALIB*21-24. Remember that stations must be analyzed in ascending 
order (the same order in which they are retrieved from the data library). 

Step 1. Fill out information sheets (see figs, 11 and 12). 

Step 2. Create data file for program use and execute program RXWTHR. 

123456 78 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

9 RUM, 

@ASG,A FIREDATALIB*PROGRAMS. 

@ASG,A FIREDATALIB*21-24. 

9USE 2., FIREDATALIB*21-24. 

@ASG,UP PHILIP. 

@USE 15., PHILIP. 

QXOT FIREDATALIB*PROGRAMS.GETDATA2. 

24020755 24020778 

24300260 24300277 
@EOF 
@PREE 2. 
@ASG,T 2. 

@ASG,A CSSG*R1LIB. 
@XQT CSSG*R1LIB. RXWTHR 

Step 3. Transfer information from user sheet(s) to machine readable formats. 
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Total \umber Different Stations in This Run ^ 

User's Same COLONY ,5ubunit GLACIER NP 

1 / 
Project- 



,Unit 



NPS 



Brief fire weather summaries for West Glacier headquarters 



Fire heather Station Information Station name I-TEST GLACICR 

T / 

46 



\ 0> 240207 



Elevation 3200 
Fuel model- 



f t . .Latitude 



.Climate class 



, Slope c 



>ear begin 1955 



.Last frost -_ 

.Year end 1978 



5/ 
_,urass type - annual 

.Date begin 06QI 



, perennia 1 



,Date end 



0915 



Site Ad i us tme nt Factors I if any) : 

\spect I, l=north, 2=east, 3= south, 4=west ), Site elevation 

Canopy cover < l=open, 2 = closed ) 



ft 



cm 
cm 
cm 
cm 



D u f f / So 1 1 Ho r i -_on__| n if o r ro a 1 1 o n ( i f D u _ff__ _Mp_i s t u r e 5 e j. e c t e d . S e e t J s i e_r ' s G ui d e , j _ ap p e n U i _.x __ DJ_ : 

La_ye_r Duff /Soil Type Thickness 

If Use up to 80 characters 

2/ See User's Guide, appendix C 

_3/ See User's Guide, appendix B 

4/ For NFDRS indices only. See User's Guide, 
appendix A. 



cm 



I 5/ For NFDRS indices only 

* * * TT* ********************* , K*it 



Summary T_abl_e_(s j____Re_ques_te_d (Select up 
_ State of the weather 
j{ Temperature (degrees F} 
Relative humiditv ( a ,) 



15) . 



Jt_Wind direction (8 point} 
Wind speed (mi/h) 



Max temperature (24 h, deg. F) 
Min temperature i24 h, deg. F) 
Max relative humidity (24 h, i 



_Min relative humidity I 24 h, " 
Precip duration (last 24 h) 

J>recip amount (24 h, 0,01 in) 
1 hour fuel moisture (%) 
10 hour fuel moisture ( D o) 

_NFDRS ERG 

_NFDRS BE 
Duff Moisture (24 h average, 



Co-occurrence Table(s) Requested (If 3-way table desired leave last space blank. If 



1 


selected, Wind Direction must always be 


listed 


first. ) ; 


WIND DIRECTION w i t h WTND SPEED 


with 




*> 


TEMPERATURE with RELATIVE HUMIDITY 


with 


WIND SPEED 


3 


with 


with 




4 


wi t h 


with 




5 


with 


with 





Figure 11.--RXWTHR user informatioii sheet for RXWTHR demonstration using the West Glacier 

weather station. 



Total \umber Different Stations in This Run 
User 1 5 Same BRADSHAW .Subunit 



02 



LAB 



.Unit INT EXP STN 



Proiect 



1 ' 



Demonstration of RXWTHR output for Philipsburg Ranger Station 



Fire he-uher 5_t J_ti on _ Information Station 



name 



PHILIPSBURG RS 



Elevation 5280 ft. .Latitude 



46 .Climate class 



, No._ 

.Slope class 



243002 



Fuel model 



-I/ 



.Last t'rost 

tear begin I960 .Year end 



57 



S/ 



, perennial 



1977 



__, Grass type : annual 

.Date begin Q501 ,p a te end Q931_ 



S i tg_\d just men t Factors Uf any) 

Aspect I l=north, 2=east, 3=south, 4=west ), Site elevation_ 

Canopv cover ( l = open, 2=closed 1 



Duff/Soil Hori-on Information (if Duff Moisture selected. See User's Guide, appendix PI: 
Layejr Ou Ft'/ Soil Type Thickness 



ft 



*************** 



cm 



cm 
cm 



cm 



cm 



]_/ Use up to 80 characters 

y See User's Guide, appendix C 

3/ See User's Guide, appendix B 

/ For NFPRS indices only. See User's Guide, 
appendix A 



5/ For NFDRS Indices only 

***T* ****************#** 



rut****************-**** 



Summary Tableisl Requested (Select up to 15] : 

State of the weather 

X Temperature ^degrees F) 

Relative humidit) ( j 

__x_Wind direction (& point) 
Kind speed (mi/h) 



Max temperature (24 h, deg. F) 
Min temperature (24 h, deg, F) 
Ma.x relative humidity (24 h, ",) 



_Min relative humidity (24 h, %) 
^Precip duration (last 24 h) 
_Precip amount (24 h, 0.01 in) 
_1 hour fuel moisture ( &) 
_10 hour fuel moisture (%) 
NFDRS ERG 



NFPRS BI 



_Duff Moisture (24 h average, 



o-occurrence Table(s) Requested (If 2-way table desired leave last space blank. If 
selected, >Vind Direction must always be listed firstTTT '' ~ 



1 WIND DIRECTION 



with WIND SPEED 



with 



TEMPERATURE 



_Wlth RELATIVE HUMIDITY Wlt h KIND SPEED 
- Wlth _ With 

with 



with 



with 



Figure 12.- RXWTHR user information sheet for RXWTHR demonstration using Philipsburg weather 

station. 



02 
STATION 

WEST GLACIER 240207 3200 46 3 
YEARS 1955 1978 
DATES 0601 0915 
ACTIVITY 

BRIEF FIRE WEATHER SUMMARIES FOR WEST GLACIER HEADQUARTERS 

1 2 34567 
1234567890123456789012345678901234567890123456789012345678901234567890 

SUMMARY 

TEMPERATURE 

WIND DIRECTION 

END 

CO-OCCUR 

WIND DIRECTION WIND SPEED 

TEMPERATURE RELATIVE HUMIDITY WIND SPEED 

END 
RUN 
02 

STATION 

PHILIPSBURG RS 243002 5280 46 3 
YEARS 1960 1977 

DATES 0501 0931 

ACTIVITY 

DEMONSTRATION OF RXWTHR OUTPUT FOR PHILIPSBURG RANGER STATION 

RUN ____________, _______ 

Step 4. Terminate program execution. 

(3EOF 

@FIN 

RXWTHR Output Interpretation 

Identical output (format) would be produced for each station by the above demonstration, 
For brevity, only output from the second station, Philipsburg, will be discussed. 
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SUTLMtt 1ABLIS 

Relanve frequency of occurrence of selected . orn 

dard deviation, median value, and value ranges are tabulated foi 10 day ana y , 

iSWkvW^^^^ 

class. Thib format is displayed in figure 2. 

of these tables should be quite straightforward; they are simply summary 
al values. In both table,, cu.ul.tive frequencies W Compute d by the 
addition of relative frequency values. For example, in May (fig. 1) 85.8 percent ^ L L P e * 
cent I 16 3 percent * 15?7 percent + 12.7 percent) of the days in the sample had temperatures 
of less than 70 degrees. Similarly, only 14.2 percent of the temperature readings weic 70 
degrees or higher. 

When obtaining summaries of NFDRS indices the user should consider excluding the first 21 
days of a year's record. The reason is that at least that much time is needed for the iuol 
moisture models to predict "converging" values for the 100- and 1000-hour fuels. (John E. 
Deeming, personal communication) 

CO-OCCURRENCE TABLESTWO-WAY 

In the two-way tables, 2 months of analysis are displayed on each page of output. One 
page (May and June) is depicted in figure 3. Values within the tables are percent frequency 
of" co-occurrence of selected class values of each parameter. In May, the most frequent co- 
occurrence condition is a windspeed of 3 to 5 mi/h from the southwest (9.7 percent). Most 
wind-speed values were in the range of 3 to 8 mi/h as seen from the columnar frequency totals. 
This information should be used for general knowledge of an area's weather and to screen and 
identify periods that have high probabilities of meeting potential prescription conditions. ^ 
After a' prescription occurrence examination by RXBURN, the information could be used in making 
new prescriptions that have high probability of occurring. 

CO-QCCURRFNCL TABLES THREE-WAY 

One month of analysis is tabulated on each page of output for three-way tables. Again, 
just one page is depicted in figure 4. The last parameter entered (windspeed, in this case) 
is stratified into five value classes; the first two are stratified into 10 value classes. 
In this example, the frequency of co-occurrence of a windspeed from 6 to 11 mi/h, a temper- 
ature of 59 degrees or less, and a relative humidity of 30 to 49 percent is 9,0 percent. It 
is obtained by adding the relative frequencies of all those particular stratifications. The 
number of days in the sample is given in the bottom right portion of the table. This infor- 
mation has the same potential for use as that from the two-way tables, with the added flexi- 
bility of a third parameter. 

Using Program RXBURN 

To use RX3URN, the four steps designated for RXWTHR are required with a slight variation in 
the substeps. The steps are relisted here with only the changes detailed. The structure of 
directive blocks is the same as for RXWTHR; they must begin with the Number of Stations card and 
end with a RUN card. 

Step 1. Fill out information sheets for RXBURN (fig. 3). 

All values up to prescription factor selections are identical to those in the RXWTHR 
sheet. A discussion on defining prescriptions for this section is given in appendix E. 

Step 2. Create data file and load program RXBURN. 

This is done exactly the same as in RXWTHR except the command @XQT CSSG*R1LIB. RXWTHR is 
altered to read: 
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@XQT CSSG*R1LIB.RXBURH 

Stop 3. Transfer information from user's sheet to machine readable formats, 

Items (aj through (d) are exactly the same as in RXWTHR. The changes in this step are 
outlined and exemplified here. In RXBURN, there is no restriction on the number of months 
analysed in one run. 

e. User identification cards. --RXBURN allows the information from the second line of 
the sheet to be entered into the program for documentation. The first card of this sequence 
instructs the program that the user's name, district, and forest follow on the next card, It 
consists of the word IDENTIFY beginning in the first column of the card. The information card 
follows with the user's name in column 1 through 15, the district's name in 21 through 35, and 
the forest's name in 41 through the 55th column. 



12 3456 
123456789012345678901234567890123456789012345678901234567890 

IDENTIFY 

WILLIAM SHAW TROUT LAKE RD CASCADE NF 

f. Prescription condition cards. --The first card of this sequence instructs the program how 
many prescription cards follow. The number (nnj may range from 1 to IS; the card consists of the 
word PRESCRIBE beginning in the first column, and the value nn in columns 11 and 12, right- 
justified. This card is then followed by the prescription condition cards. These cards are 
transferred directly from the user information sheet. The parameter name begins in column 6, the 
preferable RX minimum values are entered in column 33 through 35; preferable RX maximum values in 
columns 38 through 40; acceptable RX minimum values in columns 43 through 45; and acceptable maxi- 
mum values in the 48th through SOth column. All values must be right-justified. Values for wind 
direction should be entered as eight points in a clockwise direction; for example, southwest to 
northwest, northwest to northeast; and east to west. 



12345 

12345678901234567890123456789012345678901234567890 

PRESCRIBE 04 

WIND SPEED 6 7 4 10 

TEMPERATURE 58 59 55 63 

DAYS SINCE LAST PRECIP 7749 

WIND DIRECTION NW NW W N 

g. The RUN card. --The RUN card serves the exact purpose of that in RXWTHR and must be 
the last card in each directive block. 

Step 4. Terminate program execution. 

The program is terminated the same as RXWTHR was terminated. 

In Step 3, the order in which items (b) through (f) are entered is not restricted as long 
as the order within each item is consistant with the directions stated above. This is also 
true with input for RXWTIIR. 
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Program RXBURN Error Messages 

1. No recognisable input cards - incorrect input stream - program aborted. 

Cause: Misspelled control word on a control card, or control card word not beginning in 
column 1, 

2. At least one option card is incorrect or option word started in the wrong column - 
program aborted. 

Cause. Misspelled option after 'the PRESCRIBED control card, or option card not beginning 
in column 6. 

3. Difference between site and AFFIRMS station greater than 2,000 feet, no adjustments made. 

Cause: Self-explanatory, 

Program RXBURN Output Quantities 

For each station and/or prescription analysis in a single RXBURN run, approximately five 
pages of output are produced, with a maximum of nine in a full 12-month analysis. 

Program RXBURN Cost Estimates 

For a typical RXBURN run (one station, 5 to 6 months analysis, three to four prescription 
parameters), a user can expect the costs and resource use outlined below. 

Run priority Cost CAU (sec) Resource time (sec) 

DEMAND $4.00 to $5.00 8 to 9 10 to 11 

L $2.00 to $2.00 8 to 9 10 to 11 

N $1.00 to $1.50 8 to 9 10 to 11 

P $0.75 to $1,00 8 to 9 10 to 11 

RXBURN Demonstration 

Again, combining steps 1 through 4 results in the execution of a simple run of RXBURN. In 
this demonstration, one station is used, although more could be analyzed in a manner identical to 
that shown in RXIVTHR. It is again assumed that the data resides on file FIREDATALIB*21-24. 
Remember, the weather data are read from logical unit 15. 

Step 1. Fill out RXBURN information sheet (see fig. 13). 
Step 2. Create data file and load program RXBURN. 



26 



RXBURN USER INFORMATION SHEET 



Total Number Different Stations in This Run 01 



User's Name BRADSHAW , Subunit FIRE LAB .Unit INT EXP STH 

Project-'' Demonstration of RXBURN output for Philipsburg Ranger Station 

Fire Weather Station Information' Station name PHILIPSBURG RS $ 0t 243002 

Elevation 5280 ft. , Latitude 46 .Climate class- 3 .Slope 



Fuel model- .Last frosty .Grass type^A annual , perennial 

Year begin I960 ,Year end 1977 .Date begin 0501 ..Date end npi 

Site Adjustment Factors fif any): 

Aspect ( l=north, 2=east, 3=south, 4=west ), Site elevation ft. 

Canopy cover ( l=open, 2=closed ) 

Duff/Soil Horizon Information (if Duff Moisture selected. See User's Guide, Appendix D) : 

Laver Duff /Soil Tvne Thickness 

^."/.r^. . " r ~.>...ry.:.~ ,'.,/,Jir; V'..i.7?. *V *"?.. ************-**-****-***'****** ************* + ********** 

* 

1 __cm * y Use up to 80 characters 

2 _____________ cm * I/ See User's Guide, appendix C 

* 

3 cm * _3/ See User's Guide, appendix B 

* 

4 cm * 4_/ For NFDPS indicies only. See User's Guide, 

* appendix A 

5 cm * Sf For NFDRS indices only 



PT e _sc rip t ion F act or_ Sej. e ct_ig_n s (C h ec k_a-nd s et 1 im i t s_ for^up^ t o 1 5 factors.) : 

Pr'e f er abj^e^ _Rx Accept able^j^x 
Fj/c^tor Minimum Maximum Minimum Maximum 

State of the weather 



jj_Temperature Cdeg. F) 65 75 60 80 

_x_Relative humidity (%) 20 30 20 55 

Wind direction (S point) 



_x_JVind speed (mi/h) 4 9 15 

Max temperature (24 h , deg . F ) 



Min temperature (24 h, deg. F) , . . 
_Max relative humidity (24 h, "0 . . 
Min relative humidity (24 h, %) . , 



Days since last precipitation.. 
_Precip amount (24 h, 0.01 in). 

1 hour fuel moisture ("o) 

10 hour fuel moisture (%) 



_MFDRS ERG. , 
N'FDRS BT. .. 



Duff moisture (24 h average, 



Figure 13.--RXBURN User Information Sheet for RXBURN demonstration using Philipsburg weather 

station. 



1234 
1234567890123456789012345678901234567890, 

@RUN, . . . 

@ASG,A FIREDATALIB*PROGRAMS. 

@ASG,A FIREDATALIB*21-24. 

@USE 2. ,FIREDATALIB*21-24. 

@ASG,UP PHILIP. 

@USE 15, ,PHILIP. 

@XQT FIREDATALIB*PROGRAMS.GETDATA2 

24300360 24300277 
@EOF 
@FREE 2. 
@ASG r T 2. 
@ASG,A CSSG*R1LIB. 

CSSG*R1LIB.RXBURN 



Step 3. Transfer information from user's sheet to machine readable formats 
01 

STATION 

PHILIPSBURG RS 243002 5280 46 3 

YEARS 1960 1977 

DATES 0510 1101 

ACTIVITY 

DEMONSTRATION OF RXBURN OUTPUT FOR PHILIPSBURG RANGER STATION 

IDENTIFY 

LARRY BRADSHAW FIRE LAB MONTANA 

PRESCRIBE 03 

TEMPERATURE 65 75 60 80 

RELATIVE HUMIDITY 20 30 20 55 

WIND SPEED 4 9 15 

RUN 
tep 4, Terminate program execution. 

@EOF 
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RXBURN Output Interpretation 

Output from this demonstration is displayed in figures 5 through 8. 
SUMMARY PAGE 

Figure 5 depicts the format and values of the summary page of RXBURN output. Part 1 
consists of input information: district name, user name, station information, sample size, and 
a prescription condition summary. Part 2 is composed of a prescription occurrence summary, It 
details the mean number of days per season of each prescription type, the frequency of occurrence 
for each typo, and the highest 10-day period and month frequency of each prescription type. In 
this example, 8 days per year were preferable, 51 were acceptable, and 80 or 58 percent of the 
season had unacceptable conditions, July had the highest acceptable prescription probability at 
45 percent, and the highest 10-day period begins July 1, with a probability of 56 percent, Two 
10-day periods, June and October 1-10, had equal probabilities of preferable conditions at 10 
percent, while the least likely period for a preferable or acceptable condition to occur is the 
10-day period beginning May 1, with a probability of 22 percent of (100 percent - 78 percent) . 

Figure 6 shows the seasonal progress of prescription occurrences. The sample size, mean 
number of days, total count of occurrences, and occurrence frequency are given for each prescrip- 
tion type by 10-day, monthly, and seasonal stratifications. This information should be used to 
follow the historical pattern of prescription occurrences throughout the season. It may be help- 
ful to plot the occurrence frequencies as they change over the season to obtain a better feel for 
the progression. 

Figure 7 is an example of the output from the run length summary analysis of RXBURN. The 
mean run length (number of days in a row of the same prescription type), and the values at the 
25th, 50th (median) ( and 75th percent! le levels are given. Stratification is by 10-day and 
monthly periods. This information gives an indication of the length of time that favorable or 
unfavorable conditions tend to persist. In this demonstration, the period July 1 to July 10 
averaged a preferable run length of 2 days. At the 25th percentile level, the value was 1 day, 
2 days was the median, and the value was 3 days at the 75th percentile level. This means that 
there is a 50 percent chance of a preferable run lasting at least 2 days and a 25 percent chance 
of one lasting at least 3 days. 

The last section of output from RXBURN is the probability of a prescription type occur- 
ring in 1, 2, or 3 days, given an initial prescription type. Stratification is by month. 
Three months are displayed on each page of output. An example of this output is shown in 
figure 8. In this example, for a day in May with acceptable conditions, there is a 12 percent 
chance that tommorow' s conditions will also be acceptable, a 9 percent chance that conditions 
2 days hence will be acceptable, and a 7 percent chance for acceptable conditions on the third 
day. If today's conditions, in ths example, are unacceptable, there is an 80 percent chance 
that they will be unacceptable the next day, and only a 2 percent chance of their being acceptable, 

System Options 

ANALYZING MORE THAN ONE STATION 

Up to 99 stations may be analyzed in a single execution of either RXWTHR or RXBURN by 
building station data files (stations in ascending order) and "stacking" directive blocks (one 
for each station). This is illustrated in Step 3 of the RXWTHR demonstration. Each directive 
block must begin with a Number of Stations card and end with a RUN card. Information entered 
once within a directive block (for example, prescription condition, summary table option, or 
activity information should not be reentered in following blocks unless the information changes. 

STATIONS WITH MORE THAN 5 MONTHS WEATHER DATA 

If two or more runs of RXWTHR are required to obtain a full season analysis, set the 
Number of Stations card value to one (01) and stack two (or three) directives blocks with the 
dates of analysis cards so that 5 months are analyzed in the first run. The dates of analysis 
cards in following directive blocks should have different seasonal dates so that the entire 
season is covered in two or three runs. 
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MULTIPLE PRESCRIPTION CONDITION ANALYSIS OF ONE STATION 

Primary use of this option is to analyze two or more prescription conditions for a 
single station in one program execution. It may also be used to get more than five co- 
occurrence tables in RXWTHR. 

An initial run of the program stores computed and observed parameter values; subsequent 
runs have access to stored values, resulting in computer cost and time savings. Values are 
stored by entering a "T" in Column 48 of the first directive block's station information card, 

Stored computations are used by the entry of a "T" in column 46 of the second directive 
block's station information card (leave column 48 blank). The third and following directive 
blocks do not require' than any of the station information series be entered. There must be 
one directive block for each prescription occurrence analysis. The following example shows 
how to analyze four prescription conditions at one station. 

~1 ~~2 3 ' ~4 " 5 ..... ~ ............... 7 

1234567890123456789012345678901214567890123456789012345678901234567890 

01 

STATION (Steps 1 and 2 are assumed to have been completed) 

ANY STATION 123456 2345 45 H B P 3 0615 T 
YEARS 1900 1999 
DATES 0101 1215 
ACTIVITY 

DEMONSTRATION OF MULTIPLE PRESCRIPTION ANALYSIS 
IDENTIFY 

SANDRA LEE FERN R.S. ASPEN N.F. 

PRESCRIBE 01 

NFDRS ERG 24 1 10 

RUN 
01 
STATION 

SAME STATION 123456 2345 45 H 3 P 3 0615 T 
YEARS 1900 1999 
DATES 0101 1215 
PRESCRIBE 01 



NFDRS ERC 
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RUN 

01 

PRESCRIBE 01 

WIND SPEED 
RUN 
01 
PRESCRIBE 03 

WIND DIRECTION 

WIND SPEED 

MAX TEMPERATURE 
RUN 
@EOF 
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18 



10 



20 



S W S NW 

4 8 2 12 

65 72 60 78 



12 34567 
1234567890123456789012345678901234567890123456789012345678901234567890 

There is one restriction in using this optionparameters not computed in the initial run 
cannot be used in subsequent prescription conditions, This applies only to NFDRS indices, fuel 
moisture, and duff moisture. Fuel moisture is computed when either fuel moisture or an NFDRS 
index is specified in a prescription. This also means that fuel models and site adjustments 
cannot be altered on subsequent runs. 

SITE ADJUSTMENTS 

Site adjustments factors adjust fuel moisture values (and resulting NFDRS indices) to locations 
different from those of the base weather station. This option is provided because most weather 
data is taken at open valley-bottom locations while much prescribed burning is done on higher 
mountain slopes. The adjustment factors are from research for Fire Behavior Officer training 
courses offered through the USDA Forest Service. Use is detailed on user information sheets and 
in table 1. A brief example of a station information card to adjust to a site 800 feet higher 
than a station, north aspect, closed canopy is shown below. Site elevations must be within 2,000 
feet of the base station. It is assumed in this example that steps 1 and 2 have been completed. 



1234567 
1234567890123456789012345678901234567890123456789012345678901234567090 



ANY STATION 



123456 5000 45 3 3 



T 1 5800 2 
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APPENDIX A - NFDRS FUEL MODELS 

Fuel Model A 

This fuel model represents western grasslands vegetated by annual grasses and forbs. 
Brush or trees may be present, but are very sparse, occupying less than one-third of the area. 
txamplos of types where fuel model A should be used are cheatgrass and medusahead. Open 
pinyon-juiupcr, sagebrush-grass, and desert shrub associations may be appropriately assigned 
this fuel model if the woody plants meet the density criteria. The quantity and continuity 
of the ground fuels vary greatly with rainfall from year to year. 

Fuel Model B 

Mature, dense fields of brush 6 feet (1-8 m) or more in height are represented by this 
uiel^modcl. One-fourth or more of the aerial fuel in such stands is dead. Foliage burns 
readily. Model B fuels are potentially dangerous, fostering intense, fast -spreading fires. 
nils model is for California mixed chaparral, generally 30 years or older. The F model is 
move appropriate for pure chamise stands. The B model may also be used for the New Jersey 
pine barrens. 

Fuel Model C 

Open pine stands typify model C fuels. Perennial grasses and forbs are the primary ground 
Fuel, but there is enough needle litter and branchwood present to contribute significantly 
to the Fuel loading. Some brush and shrubs may be present, but they are of little consequence. 
Situations covered by fuel model C are open, longleaf, slash, ponderosa, Jeffrey, and sugar 
p-ine .stands. Some pinyon- juniper stands may qualify. 

Fuel Model D 

This fuel model is specifically for the palmetto-gallberry understory-pine overstory 
association of the southeast coastal plains. It can also be used for the so-called "low 
pocosins" where fuel model might be too severe. This model should only be used in the south- 
cast because of a high moisture of extinction characteristic. 

Fuel Model E 

Use this model after leaf fall for hardwood and mixed hardwood -conifer types where the 
hardwoods dominante. The fuel is primarily hardwood leaf litter. The oak-hickory types are 
best represented by fuel model B, but E is an acceptable choice for northern hardwoods and 
mixed forests of the southeast. In high winds, the fire danger may be underrated because 
rolling and blowing leaves are not accounted for. In the summer, after the trees have leafed 
out, fuel model E should be replaced by fuel model R. 

Fuel Model F 

Of the 1972 NFDRS fuel models, only the application of fuel model F has changed. Model 
F now represents mature closed chamise and oak brush fields of Arizona, Utah, and Colorado, It 
also applies to young, closed stands and to mature, open stands of California mixed chaparral. 
Open stands of pinyon-juniper are represented. Fire activity, however, will be overrated at 
low windspeeds and where ground fuels are sparse. 

Fuel Model G 

Fuel model G is used for dense conifer stands where there is a heavy accumulation of 
litter and downed woody material. Such stands are typically overmature and may also be suffering 
insect, disease, wind, or ice damage- -natural events that create a very heavy buildup of dead 
material on the forest floor. The duff and litter are deep and much of the woody material 
Ls mor-e than 3 inches (7.62 cm) in diameter. The undergrowth is variable, but shrubs are 
usually restricted to openings. Types meant to be represented by fuel model G are hemlock - 
Sitka spruce, coast Douglas-fir, and windthrown or bug-killed stands of lodgepole pine and spruce, 
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Fuel Model H 

The short-needled conifers (white pines, spruces, larches, and firs) are represented by 
fuel model H. In contrast to model G fuels, fuel model H describes a healthy stand with sparse 
undergrowth and a thin layer of ground fuels. Fires in H fuels arc typically slow spreading 
and are dangerous only in scattered areas where the downed woody material is concentrated. 

Fuel Model I 

Fuel model I was designed for clearcut conifer slash where the total loading of materials 
less than 6 inches (15.24 cm) in diameter exceeds 25 tons/acre (5.6 kg/m 2 ) . After settling 
and the fines (needles and twigs) fall from the branches, fuel model I will overrate the fire 
potential. For lighter loadings of clearcut conifer slash, use fuel model J, and for light 
thinnings and partial cuts where the slash is scattered under a residual overstory, use fuel 
model K, 

Fuel Model J 

This model complements fuel model I. It is for clearcuts and heavily thinned conifer 
stands where the total loading of materials less than 6 inches (15.24 cm) in diameter is less 
than 25 tons/ acre (5.6 kg/m 2 ). Again, as the slash ages, the fire potential will be overrated. 

Fuel Model K 

Slash fuels from light thinnings and partial cuts in conifer stands are represented by 
fuel^model K. Typically the slash is scattered about under an open overstory. This model 
applies to hardwood slash and to southern pine clearcuts where the loading of all fuel is less 
than 15 tons/acre (3.36 kg/m 2 ). 

Fuel Model L 

This fuel model is meant to represent western grasslands vegetated by perennial grasses. 
The principal species are coarser and the loadings heavier than those in model A fuels. Other- 
wise, the situations are very similar; shrubs and trees occupy less than one-third of the 
area. The quantity of fuel in these areas is more stable from year to year. In sagebrush areas, 
fuel model T may be more appropriate. 

Fuel Model N 

This model was constructed specifically for the sawgrass prairies of south Florida. It 
may be useful in other marsh situations where the fuel is coarse and reedlike. This model 
assumes that one-third of the aerial portion of the plants is dead. Fast-spreading, intense 
fires can occur even over standing water, 

Fuel Model 

The model applies to dense, brushlike fuels of the southeast. fuels, except for a 
deep litter layer, are almost entirely living in contrast to B fuels. The foliage burns 
readily except during the active growing season. The plants are typically over 6 feet fl 8 m) 
tall and are often found under an open stand of pine. The high pocosins of the Virginia and 
the North and South Carolina coasts are the ideal of fuel model 0. If the plants do not meet 
the 6-foot (1,8-m) criteria in those areas, fuel model D should be used. 

Fuel Model P 

Closed thrifty stands of long-needled southern pines are characteristic of P fuels A 
i? A + K t0 i- c ^ la y er o U ? htl y compacted needle litter is the primary fuel.' Some 
small diameter branchwood is present, but the density of the canopy precludes more than a 
scattering of shrubs and grass. Fuel model P has the high moisture of extinction characteristic 
of the southeast. The corresponding model for other long-needled pines is U. cna * a enstic 
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Fuel Model Q 

Upland Alaskan black spruce is represented by fuel model Q. The stands are dense, but have 
frequent openings filled with usually inflammable shrub species. The forest floor is a deep 
layer of moss and lichens, but there is someneedle litter and small-diameter branchwood. The 
branches are persistent on the trees, and ground tires easily climb to the tree crowns. This 
fuel model may be useful for jack pine stands in the Lake States. Ground fires are typically 
slow spreading, but a dangerous crowning potential exists. Users should be alert to such 
events and note those levels of spread component and burning index when crowning occurs. 

Fuel Model R 

This fuel model represents the hardwood areas after the canopies leaf out in the spring. 
It is provided as the off-season substitute for E. It should be used during the summer in all 
hardwood and mixed conifer-hardwood stands where more than half of the overstory is deciduous. 

Fuel Model S 

Alaskan or alpine tundra on relatively well-drained sites is the S fuel. Grass and ^ low 
shrubs are often present, but the principal fuel is a deep layer of lichens and moss. Fires 
in these fuels are not fast spreading or intense, but they are difficult to extinguish. 

Fuel Model T 

The bothersome sagebrush-grass types of the Great Basin and the Intermountain West are 
characteristic of T fuels. The shrubs burn easily and are not dense enough to shade out 
grass and other herbaceous plants. The shrubs must occupy at least one-third of the site or 
the A or L fuel models should be used. Fuel model T might be used for immature scrub oak 
and desert shrub associations in the West, and the scrub oak -wire grass type in the Southeast. 

Fuel Model U 

Closed stands of western long-needled pines are covered by this model. The ground fuels 
are primarily litter and small branchwood. Grass and shrubs are precluded by the dense canopy, 
but occur in the occasional natural opening. Fuel model U should be used for ponderosa, 
Jeffrey, sugar, and red pine stands of the Lake States. Fuel model P is the corresponding 
model for southern pine plantations. 
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APPENDIX B - NFDRS SLOPE CLASS DEFINITIONS 

Slope class Slope (.percent) 

1 0-25 

2 26-40 

3 41-55 

4 56-75 

5 75 and above 
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APPENDIX C - NFDRS CLIMATE CLASS DEFINITIONS 



NFDRS 

climate 

class 



Thornthwaite 
humidity 
province 

Arid 



Semiarid 



1 



Subhumid 
(rainfall 
deficient in 
summer) 



Subhumid 
(rainfall 
adequate in 
all seasons) 

Humid 



Wet 



Characteristic 



Desert (sparse 
grass and scat- 
tered shrubs) 



Steppe (short 

grass and shrubs) 



Savanna (grass- 
lands, dense 
brush and open 
conifer forests) 



Savanna (grass- 
lands and open 
hardwood forests) 



Forests 



Rain forest 
(redwoods, and 

spruce-cedar- 
hemlock) 



Regions 

Sonoran deserts of west Texas, 
New Mexico, southwest Nevada, 
and western Utah; and the Mojave 
Desert of California. 

The short grass prairies of the 
Great Plains; the sagebrush 
steppes and pinyon/ juniper wood- 
lands of Wyoming, Montana, Idaho, 
Colorado, Utah, Arizona, Washing- 
ton, and Oregon; and the grass 
steppes of: the central valley of 
California. 

The Alaskan interior, the chap- 
arral of Colorado, Arizona, New 
Mexico, the Sierra Nevada foot- 
hills, and southern California; 
ponderosa pine woodlands of the 
West; and mountain valleys (or 
parks) of the northern and 
central Rockies. 

Bluesteni prairies and bluestem- 
oak-hickory savanna of Iowa, 
Missouri, and Illinois. 



Almost the entire eastern United 
States; and those higher eleva- 
tions in the West that support 
dense forests, 

Coast of northern California, 
Oregon, Washington, and southeast 

Alaska. 



Thornthwaite 1931. 
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NFDRS CLIMATE CLASS 



1 P Arid and semiarid 

2D Subhumid 

3P Subhumid and humid 

4^ Wet 



Appendix C Map. NFDRS climate class (Deeming and others 1977). 



APPENDIX D - DUFF/SOIL HORIZON STRUCTURE 

Programs RXWTHR and RXBURN offer a newly applied duff moisture model. It must be stressed 
that the model is theoretical and has not been validated by field tests. It is based upon 
theoretical properties of water vapor transport through porous mediums. The model was devel- 
oped by Michael Fosberg (Fosberg 1975) and adapted for computer usage on archival fire weather 
data by a consulting firm. Jt is offered strictly on an experimental basis in the present 
versions of RXWTHR and RXBURN. Further work and field testing are planned and resulting 
changes and calibrations will be incorporated into updated versions of the programs. 

To employ this option of programs RXWTHR and RXBURN, the user needs to provide information 
about the structure of forest floor. The floor must be broken into layers of homogeneous 
materials: duff, humus, and a boundary or stable layer. At the surface, the boundary layer 
is the air, at the bottom, the layer is considered to be a sand. The user must specify the 
type of material of each layer and the thickness of each layer, The types of materials that 
the model will presently accept and the appropriate codes to be entered in the duff information 
cards are as follows: 

Duff/ soi l_type Code 

Lodgepole pine duff 1 

Lodgepole pine humus 2 

Ponderosa pine duff 3 

Ponderosa pine humus 4 

Douglas-fir duff 5 

Douglas-fir humus 6 

Sand 7 

The thickness of each layer is to entered by the user in centimeters. At this time the 
bottom boundary layer is assumed to be sand with a thickness from 100 to 300 cm (approximately 
40 to 120 inches). This information still must be entered. The type entered for each layer 
is then used to calculate layer timelags based on physical properties of each of the layers. 
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APPENDIX E - DEFINING PRESCRIPTION CONDITIONS 

UXBURN offers the unique feature of allowing two value ranges to be entered for each 
presc-ut on factor in a prescription. They are defined as "preferable'' and "acceptable" 
b n "condUions. Preferable conditions are those ideal for meeting burning object ve 
Acceptable conditions are broader (and more often occurring) ranges of values that il still 
meet burning objectives, though not as surely as burning in preferable conditions. All value 
outside the range of acceptable conditions are considered unacceptable. 



As an example, a manager might plan to burn with a prescription that included a tempci- 
ature of 58 or 59 F., last precipitation 7 days ago, and wind out of the northwest at 
6 or 7 miles per hour. Realizing this tight prescription may have a limited chance of occur- 
ring, thev consult climatological summaries from RXWTHR and decide on an acceptable prescrip- 
tion that' will probably still meet burning objectives, and yet have a reasonable probability oi 
occurring The acceptable prescription might than include a temperature of 55 to 63 ! , , last 
precipitation 4 to 9 days ago, and a windspeed of 4 to 10 miles per hour from the west to north, 
In terms of input to RXBURN, this prescription would look like this, 

Preferrable RX Acceptable RX 



Factor Minimum Maximum Minimum Maximum 

Temperature, F. 58 59 55 63 

Days since last precipitation 77 49 

Wind speed (mi/10 67 4 10 

Wind direction NW NW W N 

Any values outside the acceptable ranges would constitute an unacceptable burning condition. 
Note that wind direction limits must be entered in a clockwise manner. 

A word of caution. Windspeed and direction can be quite unreliable except Cor the site 
at which they were collected. Consequently, the manager must rely to a large extent on 
experience and knowledge of the effects of local topography on wind in order to interpret 
RXBURN values for these parameters. 
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APPENDIX F - RXWTHR AND RXBURN PROGRAM DOCUMENTATION 

This appendix contains documentation of all routines and functions o programs RXWT1IR 
and RXBURN. It is intended to serve as a reference for programers and users interested in the 
logic and structure of the programs. It provides variable identification of values stored 
in COMMON BLOCKS. 

A general overview of the structure of RXWTHR precedes its documentation. This is followed 
by a general overview of program RXBURN and documentation of the details of RXBURN that differ 
from those in RXWT1IR. 

General Overview of RXWTHR 

Program RXWTHR computes summary tables and co-occurrence tables of user-selected fire 
weather parameters. All user input is read from logical unit S in subroutine SETUP. All 
weather data from the National Fire Weather Data Library is read in subroutine CMPUTE from logical 
unit 15. Based on user options, the program will calculate either summary tables or co- 
occurrence tables, or both. In each run of RXWTHR, a maximum of 15 summary tables and five 
co-occurrence tables can be computed. 

The general flow of the program is as follows; 

Step 1: Read number of stations to be analyzed in mainline control program RXWTHR. 

Step 2; RXWTHR calls SETUP, which reads user control sequence. SETUP then calls INTERP to 
interpret user input and set program control variables and computation levels. 
INTERP then returns to SETUP, which then calculates adjustment codes if needed 
and returns to RXWTHR. 

Step 3; RXWTltR then calls CMPUTE, CMPUTE sets initial values and begins reading weather 
data and computing needed values as per user specification and summarizes data 
by calling SUM. 

Step 4: After all data have been read and values accumulated in SUM, CMPUTE calculates 
final statistics and outputs tables through two subroutines. R1TE1 outputs the 
summary tables, RITE2 outputs the co-occurrence tables, 

Step 5: After first station has been analyzed and queued to output, CMPUTE returns to 
RXWTHR to check for another station to analyze. If there are none the 
program terminates normally. If there is another station, general flow 
starts at Step 2 again until all stations have been analyzed. 

In the following pages, the COMMON BLOCKS used by RXWTHR are first defined. Mainline 
program RXWTHR is then documented. All other routines and functions are then documented in 
alphabetical order. 

RXWTHR Common Blocks 

Program RXWT1IR uses nine labeled COMMON BLOCKS and no unlabeled COMMON BLOCKS. These are 

defined as follows: 

COMMON /FTRG/ : Contains parameters used in system fire model routine. 

VALUES SET : Values in /FIRE/ are set in subroutine FUELS and subroutine CURING, all 

variables are REAL. 
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DESCRIPTION 



Weight of dead and do 1-hour fuels (pounds per square foot) 

Weight of dead and down 10-hour fuels (pounds per square foot) 

Weight of dead and down 100-hour fuels (pounds per square foot) 

Weight of dead and down I, 000-hour fuels (pounds per square foot) 

Weight of dead and dawn 1-hour fuels + cured green fuels 

Indicates first call to subroutine ERCSC 

Depth of fuel bed (inches) 

Height of dead fuels (inches) 

Height of live fuels (inches) 

Moisture content of woody- fuels (percent) 

Moisture content of herbs (percent) 

Surface area to volume ratio herbaceous fuels 

Surface area to volume ratio woody fuels 

Surface area to volume ratio 100-hour dead and down 

Surface area to volume ratio 1-hour dead and down 

Surface area to volume ratio 10-hour dead and down 

Surface area to volume ratio 1,000-hour dead and down 

Weight of herbaceous material 

Weight of cured herbacious material 

Wind factor 

Weight of woody fuels 

Moisture of extinction of dead and down fuels (percent) 

Moisture content of cured herbaceous material (percent) 



WON /LABELS/ ; 



Contains alphanumeric labels for output. All variables are integer and 
set in BLOCK DATA. 
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AY 



: DESCRIPTION 



.1(4,17) 
-2(9) 

,3(10) 
,4(17) 

L5(17) 



Header labels for all parameters 

Header labels for wind direction 

Table labels for summary tables 

Units labels for all parameters 

Short header labels for co-occurrence tables 

Month labels (JAN--DEC) 



4MON /UM/ : Contains class level limits for all parameters. All values in /L1M/ are 
integer and are set in BLOCK DATA except for LOW, which is set in the out- 
put routines. 



RAY 



MVI 



DESCRIPTION 



'(9,17) : Contains class value limits for all parameters except STATE OF THE WEATHER 



MSW(IO) 
)W(9) 



: Contains class value limits for STATE OF THE WEATHER 
: Contains lower class limits for output tables 



JMMON /OBSERV/ : Contains values from weather observation, computed moisture values, and 
other values needed for computation 



WAYS 

DSH22) 

OB CD 

QB(2) 



OB(4) 
OB(5) 

OB (6) 



DESCRIPTION 

Contains weather observation, plus computed values 

Station number 

Year 

Month 

Day 

State of the weather 

Temperature 



43 



IOB(7) 

IOB(8) 

IOB(9) 

IOB(1Q) 

IGB(ll) 

IOB(12) 

IOB(13) 

IOB(14) 

IOB(15) 

IOB(16) 

IOB(17) 

IOB(18) 

IOB(19) 

IOB(20) 

IOB(21) 

IOB(22) 



Relative humidity 

Duff moisture (computed) 

Wind direction 

Windspeed 

Stick moisture [observed) 

Maximum temperature 

Minimum temperature 

Maximum relative humidity 

Minimum relative humidity 

Precipitation duration 

Precipitation amount 

Moisture index (1, 2, or 3) 

1-hour fuel moisture (computed) 

10-hour fuel moisture (computed or observed) 

1978 NFDRS ERG (computed) 

1978 NFDRS BI (computed) 



STRUG (5, 9) 

STRUG (1,1) 
STRUG (I, 2) 
STRUG (1 , 3) 
STRUC(I,4) 
STRUG (I, 5) 
STRUC(I,6) 
STRUC(I,7) 
STRUC(I,8) 



Contains duff/soil horizon structure information as input from SETUP. There 
are up to five horizons. 

Soil or duff type code 
Thickness of layer (centimeters) 
Particle temperature timelag (seconds) 
Particle moisture timelag (seconds) 
Bulk density layer (grams per cubic centimeter) 
Particle density (grams per cubic centimeter) 
Hydraulic conductivity (centimeters per second) 
Moisture coefficient (dimensionless) 
Temperature coefficient (dimensionless) 



IELV 
I DUFF 



Adjustment elevation code (1--3) 

Duff moisture option specification (0 - no, 1 = yes) 
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IDAYB 

ID AYE 

iSLO 

IYEARB 

IYEARE 

ELEV 

LAT 

MONTHS 

MONTI IE 

RES' EAT 

SAVE 

NUMSTA 

MXT 

MNT 

MXRII 

MNRH 

I AFT 



Day to begin sample inclusion 
Day to end sample inclusion 
Adjustment slope code (1 or 2) 
Year to begin sample 
Year to end sample 
Station elevation 
Station latitude (12) 
Month to begin sample 
Month to end sample 

Logical (LI) variable that specifies if run is to read from disk file 
created by SAVE option 

Logical (LI) variable that specifies if disk file is to be created for 
subsequent REPEAT runs 

6-digit station identifier code 
Maximum temperature - ALPHA format 
Minimum temperature - ALPHA format 
Maximum relative humidity - ALPHA format 
Minimum relative humidity - ALPHA format 
Precipitation amount - ALPHA format 



COMMON /STATS/ 



Contains statistical summary values for summary and co-occurrence tables. 
Arrays for use in summary tables are stratified by month (5); 10-day 
period of month (3); and summary table number (IS). Arrays for co-occurrence 
are stratified by table number (5); month (5); first order class number (1Q} ; 
second order class number (10); and third order class number (5) transferred 
into a vector. 



ARRAY 

MEANCS.S.IS) 

30(5,3,15) 

IXTREMC5, 3,15) 
IXYEAR(5,3,15) 



DESCRIPTION 

Accumulates values stratified as described above for summary tables. 

Accumulates sums of squares and then computed values of standard deviations 

are entered, stratified as described above. 

Holds maximum values of parameters for use in range display in summary tables, 
Holds minimum values of parameters for use in range display in summary tables 
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MEDIAN(225Q) : Holds values ordered by class for computation of median values for display 
of summary tables. Values are stratified by period of month (3), month (5), 
summary table number (151, and class interval number (12), tranformed into a 
vector. 

CUMCO(12500) : Cumulator array for co-occurrence tables, stratified as described above. 
COUNT C5, 3) : Holds sample size stratified by month and period of month, 

Program RXWTHR 

Narrative: Program RXWTHR is the mainline control program for the computation of clima- 
tological summary co-occurrence tables. It keeps track of the number of stations to be analyzed 
and the number that have been analyzed. It resets the program or terminates it according to the 
number of stations requested. 

The logic is; 

Step 1: Set run number to zero, 

Step 2; Read number of stations to be analyzed. 

Step 3: If at EOF, terminate run normally; if not, continue to Step 4. 

Step 4: Call SETUP 

Step 5: Call CMPUTE 

Step 6: Rewind disk file created by SAVE option (if used) and add 1 to the run number. 

Step 7; Check to see if last station has been analyzed. If not, return to Step 2. If so, 
rewind weather data file, set run number to zero again, and return to Step 2 to see 
if stations are to be reanalyzed with perhaps different fuel models, slope classes, 
or co-occurrence tables. 

COMMON BLOCKS USED: NONE 

A listing of all SUBROUTINES and FUNCTIONS that compose the body of RXWTHR now follows. 
They are listed in alphabetical order and follow the same format as the mainline program docu- 
mentation of RXWTHR, 

RXWTHR Ructions and Subroutines 

SUBROUTINE ADSORP(I) 

Narrative; Subroutine ADSORP is used in duff moisture calculations during periods that 
sxperienced rainfall, ADSORP is called by PERCOL, which was called by DUFF, which was called 
Dy CMPUTE. The parameter I is passed into ADSORP and is the current layer in the duff horizon 
bhat is being simulated. 



:OMMON BLOCKS used: /BLOCK2/, /OBSERV/ 
""he logic is; 

: tep 1: Solve for effective percolation rate PPP, and then set QQQ and RRR from it 
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Step 2; Call ROOT- -ROOT returns value for THETA after solving a third degree polynomial equation 

Step 3: Compute particle moisture content adsorption rate (percent moisture content per second). 

Step 4: Compute actual moisture content change (amount = rate x time). 

Step 5; Adjust THETA (volumetric content of voids) by amount adsorbed by particles. 

Step 6: Calculate effective hydraulic conductivity due to water adsorption. 

Step 7: Return to PERCOL. 

SUBROUTINE ALPHA(RF,AF) 

Narrative: Subroutine ALPHA is passed the array RF, an array containing real values. It 
ENCODES the values and if the value is zero, sets it equal to an alpha blank in an A4 format. 
This makes all output of relative frequency values alpha, for easier readability. AF is the 
array that is passed back for output. ALPHA is called by RITE1 and uses no COMMON BLOCKS. 

BLOCK DATA 

Narrative: Routine BLOCK DATA sets data constants for all LABELED COMMON BLOCKS that are 
given values through the use of DATA statements, 

COMMON BLOCKS USED: /OPSET/, /LABELS/, /LIM/ 

FUNCTION CALV(F) 

Narrative: This real function converts temperature in degrees F to temperature in degrees 
Kelvin and returns to calling program. CALV is a function from the NFDRS. 

SUBROUTINE CMPUTE 

Narrative: This routine is the major computation section of the program. It reads the 
weather data and calls all other subroutines in the process, It computes statistical summaries 
and directs output procedures by calling output routines at appropriate places in the flow. 

COMMON BLOCKS USED: /SET/, /OBSERV/, /STATS/, /OPSET/, /LABELS/, /FIRE/, /LIM/ 

The logic is: 

Step 1: Initialize all constants, including site adjustment factors (IADJUS), zero out all 
values in COMMON BLOCK /STATS/. 

Step 2: If fuel model is needed as per computation level (L) , match the fuel model that was 

requested (MOD) with the test array of (MODI) and get fuel bed parameters by calling 

FUEL(II). Then set the day of last frost (LFROST) to a month and day value for use 
in subroutine CURING. 

Step 3: Compute fire weather station pressure, adding 200 feet to the station elevation if 
Alaskan. 

Step 4: Begin DO LOOP of reading weather data from logical unit 15. 
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Step 5: If REPEAT is run, skip calculations and read values from disk file, Otherwise, make 
checks for correct station and to see if the day just read is within the time period 
requested. If maximum/minimum values are requested for analysis, but missing from 
data record, read another record (Step 4) , 

Step 6: Determine if a new season has started. A new season is said to have started if the year 
changes or there is a 30-day break in the records. If new season has started, set the 
logical variable INIT to .TRUE. This allows for initial values to be set. 

Step 7: Call processing option setting routine OPTFIX. 

Step 8; Call data testing routine TESTWX. If there is an error in the data, write a message 
and read the next record. 

Step 9: If DUFF MOISTURE is requested, CALL DUFF(INIT, IFIRST) 

Step 10: If no fuel moistures or NFDRS indices are requested, go to Step 20. 

Step 11; Compute solar declination by calling function DECL and passing month and day. 

Step 12: Get 1- and 10-hour fuel moistures from subroutine ONETEN. 

Step 13: Make site adjustments if requested and add to values of 1- and 10-hour fuel moistures. 

Step 14: If no NFDRS indices are requested go to Step 19, 

Step 15: Compute 100- and 1,000-hour fuel moistures by use of routines M100 or M100A, depending 
on value returned from OPTFIX. 

Step 16: Call subroutine CURING to determine state of green vegetation, 

Step 17: Call subroutine ERCSC to determine the day's NFDRS indices. 

Step 18: Assign computed values of ERC and BI to elements 21 and 22 of IOB, respectively. 

Step 19: Assign cori-ect values of 1- and 10-hour fuel moisture to elements 19 and 20 of IOB, 
respectively. 

Step 20: If computations are to be saved for next run use, write them to the disk file that is 
called 'tape 2' . 

Step 21: Call subroutine SUM, This routine does the intermediate statistical summaries.' 
Step 22: Read a new record (Step 4) , 

At this point all records for the station have been read and summed in subroutine SUM and are ready 
for final computations and output. 

Step 30: If no summary tables have been requested, skip to co-occurrence table output section. 
Otherwise, proceed. 

Step 31: Begin DO LOOP for each summary table requested (1 to NMEANS) . 

Step 32: Get summary parameter type, order, and write out header by call to RITE1 . If fuel 

model was used, write out fuel model letter and slope class. Then write out table-tops, 

Step 33: Begin DO LOOP for each month of data (1 to NMONTH) , 
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Step 34: Begin DO LOOP for each period of the month (1 to 3). 

Step 35: Set cumulative frequencies, relative frequencies, extreme values, and if cumulative 

frequency has reached 50% get median value (MED) by calling subroutine MIDPT2 (this is 
for first class interval). 

Step 36: Begin DO LOOP for each class interval beginnin at 2. (Number 1 was set in Step 35). 

Step 37; Set cumulative frequencies, etc., as in 35 till median value is determined. Then end 
Step 36 LOOP. 

Step 58: Compute period mean, standard deviation, and output period mean, standard deviation, 
median, range values, and relative frequency distribution by calling RITE1. 

Step 39: Continue LOOP from Step 34. 
Step 40: Continue LOOP from Step 33. 

Step 41: Repeat process, this time by monthly stratifications by summing three periods for each 
month, and outputting by call to RITE1. 

Step 42; Repeat entire process for each summary table requested [Step 31) 

At this point all summary tables have been computed and written to the p~int queue. Flow continues 
in computation and output of co-occurrence tables. 

Step 50: If no co-occurrence tables were requested, return to RX.WTHR. 

Step SI: Begin DO LOOP for number of co-occurrence tables requested. 

Step 52: Set control points based upon options and order of options. 

Stop S3: If this table is a 2- way table, skip to Step 70. Otherwise, proceed. 

Stop 54: Begin DO LOOP for number of months 

Step 55; Sum total number of days in all periods of the month and write out header, fuel model, 
and slope class if needed. 

Step 56; Zero out row and column totals and begin computation loops. 
Stop 57: First write out table-tops by calling subroutine RITE2. 

Step 58: Begin DO LOOP for number of rows (class intervals) for first parameter in the table 
(1 to L4 [Option 1,1]). 

Step 59: Compute an entire row of percent frequency of co-occurrence for the class values and 
output the entire row by calling RITE2. 

Step GO: Continue Step 58 LOOP. 

Step 61: Write out column frequency totals by calling RITE2 (there are two sets). 

Step 62: Compute total frequency (row and column totals should be 100%) and write totals out 
along with the total count for the month (number of days in the sample). 

Step 63: Repeat for each month in data set (Step 54). 
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Step 64: Repeat for each table (Step 51). 

The next block is for computation and output of 2-way tables of co-occurrence. 

Step 70: Begin DO LOOP for number of months. 

Step 71: Write headers, fuel model, and slope class if needed. 

Step 72: Zero out row and column arrays (ROWTOT and COLTOT) . 

Step 73: Write out table-tops. 

Step 74; Compute simple 2-way co-occurrence frequencies and output results one row at a time. 

Step 75: Compute row and column total, sample size, and write it out. 

Step 76: Repeat for each month (Step 70). 

Step 77: Repeat for each table (Step 51). 

At this point all summary and co-occurrence tables havu been computed and queued to the printer. 

Step 80: If computations were saved on disk, put end of Tile mark (1-01 ; ) at the end of the data. 
Step 81: Return to RXWT1IR, 
SUBROUTINE CURING(INIT) 

Narrative; Subroutine CURING Is designed to model the progression of green vegetation 
through the growing season from dead, through grecnup, through transition, and again to the 
dead or cured state, The subroutine was developed by Lhe Nl'URS work unit, and taken in entirety 
from the NI'DRS. Full documentation of this subroutine is available through appropriate NFDRS 
publications. 

COMMON BLOCKS USED: /OBSERV/, /I-'IRK/ 

FUNCTION DBCL (MONTH, IDAY) 

Narrative: This function computes the solar declination angle and is taken from NPHKS 
research. It is used for calculations of 100-hour fuel moisture. 

FUNCTION DLIGtrrCLAT.DECL) 

Narrative: This function computes the total hours of daylight for the station latitude 
and the solar declination angle as computed in HHCL. It Ls also used in calculations of 100- 
hour fuel moisture. OUGHT was derived from NFDRS algorithms. 

SUBROUTINE DUFF(INIT,IFIRST) 

Narrative: This subroutine controls Lhe calculation of dut'f moisture. It Js still in the 
developmental stages and should be used with caution, "The model was developed from numerical and 
analytical solutions of the diffusion forms of the mass continuity equation and the first law of 
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thermodynamics. Analytical solutions provided a fundamental framework to evaluate nonlinear 
interactions obtained in the numerical solutions. Dimensional analysis was used to define the 
relationship between the soil properties used in the model." (Fosberg 1975.) 

COMMON BLOCKS USED: /OBSERV/, /BLOCK1/, /BLOCK2/ 

The logic is; 

Step 1: Set subroutine constants. 

Step 2: If it is the first time the subroutine has been called on this run, set up physical 
characteristics of the duff horizon. These factors include (by layer): 

1. Packing ratio (BETA) 

2. Porosity (PHI) 

3. Vapor Diffusivity (NU) 

4. Thermal Diffusivity (KAPPA) 

5. Layer Moisture Timelag (TAURHO) 

6. Layer Temperature Timelag (TAUTMP) 

7. Modified Exponential Exchange Functions (temperature [CHITMP]) 

8. Modified Exponential Exchange Functions (vapor [CHIRHO]) 

9. Depth of layer (DEPTH). 

Step 3: Ef it is initial season, call PROFL1 to set up initial temperature and relative humidity 
profiles of the horizons. 

Step 4: Begin simulation. 

Step 5: Get input values of temperature and relative humidity. 

Step 6; Begin 24-hour duff moisture simulation (2-hcur increments). 

Step 7: Use harmonic analysis to model temperature and relative humidity cycles for each time 
increment. 

Step 8: If TPTAMT is greater than zero, and simulation time is less than precipitation duration, 
it is still raining. GO TO Step 10. 

Stop 9: This means that it has either stopped raining or did not rain. If EMC is less than 

FIBSAT, regular solution will be used and SWITCH is .TRUE, for the layer. Otherwise, 

moisture is lost at the above saturation point rate as defined by DTHETA. This is done 
for each layer, and skips to Step 11. 

Stop 10- Call PERCOL. This sets a percolation rate of free water through the horizons ^ which is 
P adsorped to an extent by a value of DTHETA, and contributes to the overall moisture 

content of the horizon by a volumetric amount that is determined by solution of a third 
degree polynomial equation by subroutine ROOT. This value is determined for each layer, 
and what is adsorped by one layer is not available for adsorption by the next. 

Step 11: Determine vapor pressure for each layer and set top layer boundary conditions and inter- 
mediate layers to last value of layer above. 

Step 12: Modify layer temp and density by exchange functions. 

Step 13: Recalculate vapor pressure using new layer temperature (sat v.p.). 

Step 14: Compute current saturation density. 
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Step 15: If EMC greater than FIBSAT, current density = saturation density. 

Step 16: Compute vapor pressure and layer relative humidity from equation of state. 

Step 17: If SWITCE1 of layer is .TRUE, get particle moisture by calling DUFFMC which calculates 

EQMC of particles based on temperature and relative humidity. Compute average, simulate 
next time period, until 24 hours have been simulated, then return to CMPUTE with IOB(8) = 
duff moisture. 

SUBROUTINE DUFFMC(TEMP, HUMID, EQM, TYPE) 

Narrative: Subroutine DUFFMC computes duff particle equilibrium moisture contents based 
on a regression shown in the Wood Products Handbook, using predictors of temperature and 
relative humidity. It also computes EQMC's for soil or loam boundary layers. 

FUNCTION EQMC(TEMP,RH) 

Narrative: This function provides the same information as DUFFMC, but is in the form of 
a function that does not have the capability for soil and boundary layer computations. 

SUROUTINE ERCSC(ERC,BI) 

Narrative: This computational routine determines potential fire behavior from the NFDRS 
fuel model and the current fire weather. It then computes the fire danger indices based on 
the potential fire behavior. The model is based on work of Rothermel (1972) and Albini (1976) 
and modified by Deeding and others (1977). This routine is taken directly from the NFDRS 
system and again is fully documented there. 

COMMON BLOCKS USED: /FIRE/ 

SUBROUTINE ERROR(I) 

Narrative: This subroutine is called by several routines in the system and is passed 
the value I. The routine is called when an error is detected. An error message is printed 
by means of a computed GO TO statement and processing continues or the program is aborted, 
depending on the severity of the error. 

FUNCTION FREEZE (ITMN,NUMDAY) 

Narrative: This function determines if a killing frost has occurred and is called by 
CURING, A killing frost is considered to have occurred if the minimum temperature lias been 
less than 33 degrees for a total of 5 days, or if a minimum temperature of 25 degrees or less 
occurs. This logical function is taken from the NFDRS, 

SUBROUTINE FUEL (I I) 

Narrative: This routine sets the parameters of the fuel bed for analysis or modeling by 
ERCSC. The routine contains the 20 NFDRS stylized fuel models, and assigns the correct values 
to the variables in COMMON BLOCK /FIRE/. FUEL is called by CMPUTE. 

IBROUTINE INTERP(MVALUE,MOPTIN) 

Narrative: This routine interprets what was read in SETUP, sets program directives, and 
amputation levels (L) . It is passed the alphanumeric words that were read by SETUP, matches 
.hem with its own word dictionary, and sets option parameters and computation levels. If no 
latch is found, ERROR is called with the appropraite value passed. 

OMMON BLOCKS USED; /SET/, /OPSET/ 
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FUNCTION lOUT(KSET) 

Narrative: This function examines the option number and sets a value (TOUT) for use in 
COMPUTED GO TO in the output routines. 

FUNCTION IRND(X) 

Narrative: This integer function rounds a real number then performs a FORTRAN standard 
IPIX operation on it. 

INTEGER FUNCTION ISUB1 (I, J, K, L) 

Narrative: This function converts a four-way stratification into a vector for computer 
storage on a computer that is restricted to three dimensional arrays. This particular function 
is used in computing median values in the summary array MEDIAN. 

INTEGER FUNCTION ISUB2 (I, J, K, L,M) 

Narrative: This function converts a five-way stratification into a vector for computer 
storage. This particular function is used in co-occurrence stratification of array CUMCO. 

INTEGER FUNCTION JET(I) 

Narrative: This function groups precipitation and ERG values into categories of graduated 
value. 

SUBROUTINE M100CDECL.LAT, BNDRY) 

Narrative: The routine is called by CMPUTE to calculate 100-hour fuel moistures. It is 
the preferred method of 100-hour calculations. Maximum and minimum temperatures and relative 
humidities are used to calculate extreme values of equilibrium moisture contents which are then 
integrated with PLIGHT resulting in a weighted average value that is used in the 100-hour 
computations. This routine was developed by the NFDRS work unit and includes BNDRY values. 

COMMON BLOCKS USED: /OBSERV/ 

SUBROUTINE M100A(IWTFLG, LAT, DECL, BNDRY) 

Narrative: This routine developed by the NFDRS work unit is for computing values of 100- 
hour fuel moisture when maximum and minimum temperatures and relative humidity are not avail- 
able, It uses default values of moisture recovery at night based on climate classes and 
yesterday's 100-hour fuel moisture content. It then computes a boundary condition that is 
used in computations of 1,000-hour fuel moistures. 

COMMON BLOCKS USED: /OBSERV/ 
SUBROUTINE M1000(BNDRY,MC1000, INIT) 

Narrative: This NFDRS routine computes 1,000-hour fuel moisture by using a 7-day average 
boundary (BNDRY) condition as a driver function in the calculations. The value for BNDRY is 
computed in either M100 or M100A. 

SUBROUTINE MIDPT2(M,C, X,CO, J, L, LI) 

Narrative: This subroutine computes the median value of an order ranked distribution by 
methods for large sample sizes described by Nie and others (1975). M is the computed median 
value, C is the cumulative frequency, X is the last class's relative frequency, CO is the 
sample size, J is an an/off switch, and L and LI are upper and lower class value limits. 
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SUBROUTINE MODE(RF, J, ISET) 

Narrative: This subroutine selects the modal class value from a relative frequency 
distribution. RF is the distribution, J is an on/off switch, and ISET is an option delimiter 
set in FUNCTION I OUT. 

INTEGER FUNCTION NDEX(I, J, K, L,M) 

Narrative: This function- computes the vector element of a 5 -way stratification adjust- 
ment vector. The stratifications are (1) aspect, (2) elevation, (3) slope, [4) canopy cover, 
and (5) month, 

SUBROUTINE ONETEN 

Narrative: This function computes 1- and 10-hour fuel moisture values from the daily 
fire weather observation. This routine was developed by the NFDRS and has several additional 
data integrity checks for use in this program. 
The logic is: 

Step 1: Set temperature and relative humidity adjustment factors to compute fuel bed tempera- 
ture and relative humidity based on cloud cover. 

Step 2: Calculate equilibrium moisture content by FUNCTION EQMC. 

Step 3: If stick moisture is reported and in allowable limits, use it and calculate 1-hour 
fuel moisture from 10-hour value and return. 

Step 4; Use regression equation to compute 1- and 10-hour values if stick moisture is not 
reported, 

Step 5: Return to CMPUTE. 
COMMON BLOCKS USED: /OBSERV/ 

SUBROUTINE OPTFIX 

Narrative: This NFDRS routine computes which processing option will be used for fuel 
moisture computations. IOPTIN = 1 is the preferred option and IOPTIN is the least desirable 
option (due to mimimal data), 

COMMON BLOCKS USED: /OBSERV/ 

SUBROUTINE PERCOL 

Narrative: This subrouitne is called by DUFF, in the case of precipitation, to compute 
the percolation rate of free water through the duff horizons. 

COMMON BLOCKS USED: /OBSERV/, /BLOCK2/ 

The logic is: 

Step 1: Compute precipitation rate. 
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Step 2: Set percolation rate as the minimum value of the precipitation rate of the hydraulic 
conductivity of the layer. 

Step 3: Call ADSORP to compute THETA and DTHETA for the first layer. 

Step 4: Adjust EMC and DMC as affected by THETA and DELTAM set SWITCH to FALSE. 

Step 5: Repeat for lower layers. 

Step 6: Return to DUFF. 

SUBROUTINE PROFL1 

Narrative: This subroutine is called by DUFF to initialize soil/duff horizons to initial 
season values. It is called at the beginning of each new season. The bottom layer (organic) 
and the lowest horizon (nonorgamc] are assumed to be at fiber saturation. A linear gradient 
is established from the bottom layer to the surface relative humidity value, 

COMMON BLOCKS USED: /OBSERV/, /BLOCK2/ 

The logic is: 

Step 1: Set boundary layer initial conditions. 

Step 2: Determine bottom layer subscript. 

Step 3: Assume bottom layer (organic) and lowest horizon (nonorganic) are saturated. 

Step 4: Establish linear gradients of temperature and humidity between top and bottom layers. 

Step 5: Calculate vapor pressure and vapor densities of each horizon. 

Step 6: Return to DUFF. 

SUBROUTINE RELHUM 

Narrative: This adaptation of an NFDRS routine computes the relative humidity given wet 
and dry bulb temperatures. 

SUBROUTINE RITE1 

Narrative: This is the main output subroutine for summary tables. It first calls ALPHA 
to change all output to a "A" format and then goes to the correct writing sequence via a series 
of computed GO TO statements operating on variables passed in the calling routine CMPUTE, 

COMMON BLOCKS USED: /LABELS/, /SET/, /LIM/ 

SUBROUTINE RITE2 

Narrative: This is the main output subroutine for co-occurrence tables. It first calls 
ALPHA to convert REAL values to an "alpha" format. It then goes to the proper entry point via 
computed GO TO statements operating on variables passed in the calling routine CMPUTE, 

COMMON BLOCKS USED: /LABELS/, /LIM/ 
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SUBROUTINE ROOT (NC, III, JJJ, KKK, PPP,QQQ,RRR,ADFM) 

Narrative: This routine solves an Nth degree polynomial by numeric methods. The value 
returned in ADFM(l) is the volumetric moisture content the duff void space for solution to 
free water due to rain infilitration through the layers, The routine was given to the project 
by Michael Fosberg. 

SUBROUTINE SETUP 

Narrative: The call to SETUP reads user run information cards and "sets up" the run. It 
is called by RXWTHR. It calls INTERP and returns to RXWTHR. By reading an interpreting control 
card word, the program executes proper reading sequences by means of computed GO TO statements. 

COMMON BLOCKS USED: /OPSET/, /SET/, /OBSERV/ 

The logic is: 

Step 1; Read control card word. 

Step 2: Interpret word from internal word file. If RUN card, call INTERP, otherwise execute 
proper read sequence and return to step 1. 

Step 3: Compute adjustment indices if required. 
Step 4: Return to RXWTHR. 
SUBROUTINE SUM 

Narrative: This is an intermediate computational routine for accumulating parameters 
needed for means, standard deviation, range, median, and frequency distribution calculations. 

COMMON BLOCKS USED: /SET/, /STATS/, /OBSERV/, /OPSET/ 

The logic is: 

Step 1: Compute month index and period number and increase value in sample size array COUNT 
by 1. 

Step 2; Begin DO LOOP for each summary table requested. 

Step 3: Set array pointer value K, based on summary parameter. 

Step 4: Accumulate values for mean and sum of square computations. 

Step 5: Compute interval class of summary parameter, 

Step 6: Increment correct vector element of MEDIAN by 1. 

Step 7: Set high and low values for range analysis, 

Step 8: Continue summary DO LOOP until executed once for each summary parameter. 

Step 9: If any co-occurrence tables requested, begin DO LOOP from 1 to NTABLE. Otherwise 
return to CMPUTE. 

Step 10; Set class interval and index values for each variable in current table. 
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Step 11: Increment correct element of co-occurrence vector CUMCO by 1. 
Step 12: Continue DO LOOP. 
Step 13: Return to CMPUTE. 
SUBROUTINE TESTWX 

Narrative: This NFDRS routine tests each day's weather data for errors, computes relative 
humidity of need, and checks or computes precipitation duration. It sets some potentially 
undefined variables in the IOB array when used on REPEAT runs, 

COMMON BLOCKS USED: /OBSERV/ 

FUNCTION U20FNC(SGBRT,U20,WNDFC) 

Narrative: This function reduces the 20-foot observed windspeed to wind at wind flame height 
for calculations in ERCSC. It is called from ERCSC and is fuel bed dependent. 

SUBROUTINE VAPOR (T,EWS) 

Narrative: This subroutine computes saturation vapor pressure using formulas from the 
Smithsonian tables (Smithsonian Institution 1939). It is from the NFDRS and returns pressure 
in millibars. 

General Overview of RXBURN 

Program RXBURN computes local prescription condition frequencies based on user-defined. 
prescription conditions. All user input (except for the number of stations to be analyzed) 
is read in SETUP from logical unit 5. All weather data are read from logical unit 15 in sub- 
routine CMPUTE. All output (except for error messages) is done from subroutine MEAN. RXBURN 
allows up to 15 prescription factors to be entered into a prescription. 

The general flow of the program is as follows: 

Step 1: Read number of station to be analyzed in mainline control program RXBURN. 

Step 2; RXBURN calls SETUP, which reads user control sequence. SETUP then calls subroutine 

INTERP to interpret user input, set program control variables and computation levels, 
INTERP then returns to SETUP, which then computes adjustment codes if needed and returns 
to RXBURN. 

Step 3: RXBURN then calls CMPUTE. CMPUTE sets initial values and begins reading weather data 
and computing needed values as per user specification. 

Step 4: CMPUTE then calls SUMMER, which tests to see if the prescription factors just read/ 

computed constitute a preferable condiiton, an acceptable condition, or an unacceptable 
condition. SUMMER then accumulates values needed for mean values, run lengths, and 
persistence probabilities. SUMMER then returns to CMPUTE, which reads the next day's 
weather , 

Step 5: When all weather data has been read, analyzed, and summed, CMPUTE calls MEAN. Sub- 
routine MEAN computes all final statistics and writes all output to the printer queue. 
MEAN then returns to CMPUTE. 

Step 6: CMPUTE returns to RXBURN to check for another station to analyze, If there is none, 
the program terminates normally. If there is another station, general flow returns 
to step 2. This is repeated until all stations have been analyzed, 
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In the following pages, the COMMON BLOCKS used by RXBURN are first defined. Mainline 
program RXBURN is then documented, followed by all other subroutines and functions in alpha- 
betical order. 

RXBURN Common Blocks 

Program RXBURN uses seven labeled COMMON BLOCKS and no unlabeled COMMON BLOCKS. All values 
are set in the programs and none in DATA statements. The COMMON BLOCKS are defined as follows; 



COMMON /BLOCK1/ 
COMMON /BLOCK2/ 
COMMON /FIRE/ 
COMMON /OBSERV/ 
COMMON /SET/ 

VALUES SET 

ARRAYS 

ACT (20) 

DIST(4) 

FOREST (4) 

IRXOP(IS) 

LAMAT(15,2,2) 

LIMIT (15, 2,2) 
NAME (4) 
NAMSTA(4) 
NRXOP(15) 

VARIABLES 

ADJUST 

ELEV 

IASP 

ICOV 

IDAYB 



; This common block is the same as /BLOCK1/ of RXWTIIR. 

: This common block is the same as /BLOCK2/ of RXWTIIR. 

; >This common block is the same as /FIRE/ of RXWTIIR. 

: This common block is the same as /OBSERV/ of RXWTIIR. 

: Values needed for the setting up and control of Program RXBURN are held 
in this common block. It is basically a combination of values from 
/OPSET/ and /SET/ of RXWTHR, plus variables needed in RXBURN. 

: Most values in this block are set in subroutine SETUP and subroutine 
INTERP. 

: DESCRIPTION 

: Contains alphanumeric activity information (20A4 elements). 

: Contains alphanumeric name of forest district (4A4) , 

: Contains alphanumeric name of forest (4A4) ( 

: Contains pointer values for extracting values from the array IOB. 

: Contains alphanumeric prescription limits stratified by prescription param- 
eter, preferable prescription (RX) minimum limit, preferable RX maximum 
limit, acceptable RX minimum limit, and acceptable RX maximum limit. 

: Contain integer limits for prescription, stratified as in LAMAT. 

: Contains alphanumeric user name (4A4) , 

: Contains alphanumeric station name (4A4) . 

: Contains the order values stored in IRXOP. 

: DESCRIPTION 

; Logical switch for inclusion of adjustment option. 

: Elevation of base fire weather station (feet). 

: Aspect adjustment code for use when ADJUST is true, 

: Canopy cover adjustment code for use when ADJUST is true. 

: Day of year to begin inclusion of days in sample. 
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IDAYE : Day of year to end inclusion of days in sample. 

IDUFF : On/off switch that indicates if duff moisture is to be calculated, 

Program RXBURN 

Narrative; Program RXBURN is the mainline control program for the assessment of local 
prescription condition frequencies. _ It keeps track of the number of stations to be analyzed 
and the number that have been analyzed. 

The logic is: 

Step 1; Set run number to zero. 

Step 2: Read number of stations to be analyzed. 

Step 3: If at EOF, terminate the run normally. If not, continue to Step 4. 

Step 4: Call SETUP. 

Step 5: Call CMPUTE. 

Step 6: Rewind disk file created by SAVE option (if used) and add 1 to the run number. 

Step 7; Check to see if last station has been analyzed. If not, return to Step 2* If it has 
been analyzed, rewind weather data file, set run number to zero, and return to Step 2. 

SUBROUTINE ADSORP 

Narrative: Same as RXWTIIR. 
FUNCTION CALV 

Narrative: Same as RXWTHR. 
SUBROUTINE CMPUTE 

Narrative: This routine is a major weather reader and computational block. It performs 
the exact processes that make up the balance of CMPUTE in RXWTHR, except that no statistics 
are computed. In RXBURN, the statistics and output are handled in subroutines SUMMER and MEAN. 

COMMON BLOCKS USED: /SET/, /OBSERV/, /FIRE/ 

The logic is: 

Steps 1 through 22 are essentially the same as RXWTHR except for two items. First, alpha input 
LIMIT values are DECODED to integer format, and secondly, SUMMER is called instead of SUM at 
Step 21. 

Step 30: Call SUMMER to close out last observation strings, and ENCODE LIMIT back to alpha 
formats for output. 

Step 31: Call MEAN and put EOF on logical unit 2 (disk file). 
Step 32: Return to RXBURN. 
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RXBURN Functions and Subroutines 

SUBROUTINE CURING 

Narrative: Same as RXWTIIR. 
FUNCTION DECL 

Narrative: Same as RXWTHR. 
FUNCTION DLIGHT 

Narrative: Same as RXWTHR. 
SUBROUTINE DUFF 

Narrative: Same as RXWTHR. 
SUBROUTINE DUFFMC 

Narrative; Same as RXWTHR. 
FUNCTION EQMC 

Narrative: Same as RXWTHR, 
SUBROUTINE ERROR 

Narrative: Same as RXWTHR. 
FUNCTION FREEZE 

Narrative; Same as RXWTHR. 
SUBROUTINE FUEL 

Narrative: Same as RXWTHR. 
FUNCTION INDEX 

Narrative: Same as RXWTHR. 
SUBROUTINE INTERP 

Narrative: Same as RXWTHR with changes for RXBURN, 
FUNCTION IRND 

Narrative: Same as RXWTHR. 
SUBROUTINE M100 

Narrative: Same as RXWTHR. 
SUBROUTINE M100A 

Narrative: Same as RXWTHR. 
SUBROUTINE Ml 000 

Narrative: Same as RXWTHR, 
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SUBROUTINE MAX 

Narrative: This subroutine selects the maximum occurrence frequency for summary output 
of highest month and 10 days of frequency of prescription occurrence. 

SUBROUTING MID 

Narrative; This subroutine selects quartile values of run length as called by MEAN at 
the right time, and sets a switch that indicates value has been computed. 

SUBROUTINE MEAN 

Narrative: This subroutine makes the final calculations of accumulation arrays and does 
all program output. 

COMMON BLOCKS USED: /SET/, /OBSERV/, /STATS/, /STATS2/ 

The logic is: 

Step 1: Set initial values and call ZER02. 

Step 2: Set total counts for seasonal, monthly, and 10-day period values, and write out first 
page head and prescription entries. 

Step 3: Compute average number of days of each prescription type per season and compute 
associated frequencies and output values. 

Step 4: Select highest monthly frequency of each prescription type and write them out. As many 
as two highest months are allowed for cases of equal frequencies. 

Step 5: Select 10-day period of highest prescription frequency of each type and output. 
Again, up to two periods are allowed. 

Step 6; Write out monthly and 10-day summaries of seasonal progression of prescription 
condition occurrence. 

Step 7: Write out seasonal totals. 

Step 8; Calculate and output mean values and quartile values of run length by monthly and 
10-day periods. 

Step 0: Using a simple Markov model, compute and write persistence and transition probabil- 
ities. Stratification is by month. 

Step 10: Return to CMPUTE. 
FUNCTION MTU 

Narrative: This function sets a month pointer for use in output of the alpha month. 

SUBROUTINE ONCTEN 

Narrative: Same as RXWTHU. 
SUBROUTINE OPTFIX 

Narrative: Same as RXWTHR. 
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SUBROUTINE PERCOL 

Narrative; Same as RXWTHR. 
SUBROUTINE PROFL1 

Narrative: Same as RXWTHR. 
SUBROUTINE RELUM 

Narrative; Same as RXWTHR. 
SUBROUTINE ROOT 

Narrative; Same as RXWTHR. 
SUBROUTINE SET 

Narrative; Sets switch values for use in determination of maximum frequencies in MEAN. 
SUBROUTINE SETUP 

Narrative; Same as SETUP in RXWTHR with appropriate changes for RXBURN values. 
SUBROUTINE SUMMER 

Narrative: This routine fills the accumulating arrays contained in COMMON BLOCK /STATS/. 
The logic is: 

Step 1: Set initial valuesmonth number, period number, and increment sample size counter. 
Step 2: Check for beginning of season (INIT) and set appropriate value as needed. 
Step 3: Determine if day is "preferrable," "acceptable," or "unacceptable." 
Step 4; Increment mean value accumulator arrays. 
Step 5; Increment appropriate run length counters. 

Step 6: Determine if prescription "type" string was broken. If not go to Step 8. 
Step 7; Set appropriate run length counter arrays depending on what "type" of string was broken. 
Step 8: Set counting arrays for use in Markov model. 

Step 9: Set all today's values to holding value to become yesterday's value on the next day. 
Step 10: Return to CMPUTE. 
SUBROUTINE -TESTWX 

Narrative: Same as RXl^HR. 
FUNCTION U20FNC 

Narrative: Same as RXWTHR. 
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SUBROUTINE VAPOR 

Narrative: Same as RXWTHR. 
SUBROUTINE ZER01 

Narrative: Sets all values in COMMON BLOCK /STATS/ to zero for each run of RXBURN, 

COMMON BLOCKS USED: /STATS/ 

SUBROUTINE ZER02 

Narrative: Sets all values in COMMON BLOCK /STATS2/ to zero for each run of RXBURN, 
COMMON BLOCKS USED: /STATS2/ 
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The Iiuermountam Station, headquaiteied in Ojjden, 
Utah, is one of eight le^ional e.xpennient stations chained 
\vith pioviduig scientific knowledge to help icsouro 1 
manageis meet human needs and protect 1'oiest and mnj',0 
ecosystems. 

The Intermountam Station includes the States of 
Montana, Idaho, Utah, Nevada, and western Wyoming. 
About 23! million acics, or 85 peieem, of I he land ai cam the 
Station territory are classified as forest and lanj^eland. These 
lands include giasslands, deserts, shiublands, alpine aicas, 
and well-stocked foiests. They supply fibei for forest in- 
dustries; mmeials foi energy and industiial development; and 
water for domestic and industrial consumption. They also 
ptovide i cci cation opportunities tor millions of visitois each 
year. 

Field progiams and research work units of the Station 
aie maintained in: 

Boise, Idaho 

Bozeman, Montana (in cooperation with Montana 
State University) 

Logan, Utah (in coopeiujffbp with tltah State 

University) <' .. / 

g o ' > 
Missoula, Montana (in Epope&ktioriJVith the 

University of Montana) p. 1 ' L; 

t .- Ic"; 

Moscow, Idaho (in cooperation with t^b Univer- 
sity of Idaho) tl\ ^ r\ 
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Provo, Utah (in cooperation/with Hicham Young 
University) "*"' 

Reno, Nevada (in cooperation with the Umveisity 
of Nevada) 




